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Abstract 

Coronal loops are the building blocks of the X-ray bright solar corona. They owe their 
brightness to the dense confined plasma, and this review focuses on loops mostly as structures 
confining plasma. After a brief historical overview, the review is divided into two separate 
but not independent sections: the first illustrates the observational framework, the second 
reviews the theoretical knowledge. Quiescent loops and their confined plasma are considered, 
and therefore topics such as loop oscillations and flaring loops (except for non-solar ones which 
provide information on stellar loops) are not specifically addressed here. The observational 
section discusses loop classification and populations, and then describes the morphology of 
coronal loops, its relationship with the magnetic field, and the concept of loops as multi- 
stranded structures. The following part of this section is devoted to the characteristics of 
the loop plasma, and of its thermal structure in particular, according to the classification 
into hot, warm and cool loops. Then, temporal analyses of loops and the observations of 
plasma dynamics and flows are illustrated. In the modeling section starts some basics of loop 
physics are provided, supplying some fundamental scaling laws and timescales, a useful tool 
for consultation. The concept of loop modeling is introduced, and models are distinguished 
between those treating loops as monolithic and static, and those resolving loops into thin and 
dynamic strands. Then more specific discussions address modeling the loop fine structure, 
and the plasma flowing along the loops. Special attention is devoted to the question of loop 
heating, with separate discussion of wave (AC) and impulsive (DC) heating. Finally, a brief 
discussion about stellar X-ray emitting structures related to coronal loops is included and 
followed by conclusions and open questions. 
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1 Introduction 



The corona is the outer part of the solar atmosphere. Its name derives from the fact that, since it is 
extremely tenuous with respect to the lower atmosphere, it is visible in the optical band only during 
the solar eclipses as a faint crown (corona in Latin) around the black moon disk. When inspected 
through spectroscopy the corona reveals unexpected emission lines, which were first identified as 
due to a new elem ent (coronium) but which were later ascertained to be due to high excitation 
states of iron (e.g.. lGolub and Pasachol . ll997l : lGolub and Pasachofl1 . l200l[ ). It became then clear 



that the corona is made of very high temperature gas, hotter than 1 MK. Almost all the gas is fully 
ionized there and thus interacts effectively with the ambient magnetic field. It is for this reason 
that the corona appears so inhomogeneous when observed in the X-ray band, in which plasma at 
million degrees emits most of its radiation. In particular, the plasma is confined inside magnetic 
flux tubes which are anchored on both sides to the underlying photosphere. When the confined 
plasma is heated more than the surroundings, its pressure and density increase. Since the tenuous 
plasma is optically thin, the intensity of its radiation is proportional to the square of the density, 
and the tube becomes much brighter than the surrounding ones and looks like a bright closed arch: 
a coronal loop. 

When observed in the X-ray band, the bright corona appears to be made entirely by coronal 
loops that can therefore be considered as the building blocks of X-ray bright corona. This review 
specifically addresses coronal loops as bright structures confining plasma. It first provides an obser- 
vational framework that is the basis fo r the second part dealing with modeling and interpretation - 



Ther e have been se veral earlier books dBrav et al . 199ll;lGolub and Pasachomll997l;lAschwandenl . 
20f)4l )_ and reviews (IVaiana and Rosnerl. Il978l: Peres and Vaianal. Il99dl: Golubl. 1 19961: Aschwanden 



et l200ll : iReald . l2005h . in particular on coronal heating (jZirkerl . Il993l : ICargilll . Il995l : iKlimchukl . 



2006). that have in general a larger or different scope but include information about coronal loops 



Interested readers are urged to survey these other reviews in order to complement and fill in any 
gaps in topical coverage of the present paper. 
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2 Historical Keynotes 



First evidence of magnetic confinement came from rocket missions in the 196 0s. In particular, i n 



1965, arcmin angular resolution was achieved with grazing incidence optics ([Giacconi et aLl . ll965[ ) 



The data analysis led to the first density and temperature diagnostics with wide band filters, 
to derive high pressure in compact regions w ith intense bipolar magnetic fields and to propose 



the magnetic confinement ( Reidy et al. . 19681 ). The first coronal loop structures were identified 



properly after a rock et launch in 1968 which provided for the first time an image of an X-ray flare 



( Vaiana et al. . 19681 ). with a resolution of a few arcsec 



gineering (AS&E) program, [Vaiana et all (|l973f ) proposed a classification of the morphology of 
the X-ray corona as fundamentally consisting of arch-like structures connecting regions of opposite 
magnetic polarity in the photosphere. The classification was based on the loop size, on the physi- 
cal conditions of the confined plasma, on the underlying photospheric regions. They distinguished 
active regions, corona l holes, active regions interconnection, filament cav ities, bright points, and 
large-scale structures (jVaiana and Rosnerl . 1978; iPeres and Vaia na. 1990). 

The magnetic structuring of the solar corona is evident. However, the magnetic field lines 
can be traced only indirectly because direct measurements are feasible generally only low in the 
photosphere through the Zeeman effect on spectral lines. It is anyhow possible to extrapolate the 
magnetic field in a v olume. This was do ne to deri ve the magneti c field structure of a relatively 
stable active region bv lPoletto et al. (1975) using the Schmidt f)l964h method, under the assumption 
of negligible currents in the corona. This was useful also to derive magnetic field intensities 
sufficient for hot plasma confinement. Later on, even m ore reliable magnetic field topologies were 
derived assuming force- free fields (e.g.. ISakural 119811 ). i.e., with currents everywhere parallel to 
the magnetic field as it is expected in coronal loops. 

The rocket missions lacked good time coverage and the information about the evolution of 
coronal loops was only limi t ed, m ostly available from the Orbiting Solar Observatory-IV (OSO- 



IV) mission ( Krieger et al. . 1972h . This satellite had an angular resolution in the order of the 



arcmin and could not resolve individual loops. In 1973 the X-ray telescope S-054 on-board Skylab 
monitored the evolution of coronal loops for several months, taking 32,000 X-ray photographs with 
a maximum resolution of 2 arcsec and an extended dynamic range. It was possible to study the 
whole evolution of an active region, from the emergence as compact loops filled with dense plasma 
to its late spreading, a few s olar rotations l ater, as progressively longer and longer loops filled with 
less and less dense plasma (jGolub et It was confirmed that the whole X-ray bright 

corona consist s of magnetic loops , whose lifetime is typically much longer than the characteristic 



cooling times ( Rosner et al. . 1978f ). No exception is made by the coronal holes where the magnetic 



field opens radially to the interplanetary space and the plasma streams outwards with practically 
no X-ray emission. 

In the same mission coronal loops were detected also in the UV band at temperatures below 
1 M K bv Extreme UltraViolet CEUVI telescopes S-055 dReeves et all ll977Tl and S-082 f Tousev 



et al. Il977l ; iBartoe et <zZ.L Il977l ) . These loops are invisible in the X-ray band and many of them 
depart from s unspots , appear coa xial and are progressively thinner for progressively lower tem- 
perature ions (jFoukall . Il975l . 1976), suggesting nested loops hotter toward outer shells, although 
this suggestion has never been validated. The apparent scale height of the emission is larg er than 



that e xpected from a static model, but the loops appear to be steady for long times. iFoukall 
(Il976h proposed a few explanations including siphon flows and thermal instability of the plasma 
at the loop apex. New observations of such cool loops were performed several years later with the 
SOlar and Heliospheric Observatory (SoHO) mission and provided new details and confirmations 
(Section l3~3|) . 

A different target was addressed by the Solar Maximum Mission (SMM, 1980-1989, Bohlin 
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et aL ll98d : lActon et all fl980h ■ which included high-resolution spectrometers in several X-ray lines 
i.e. the Bent Crystal Spectrometer (BCS) and the Flat Cryst al Spectrometer (FCS), mostly devoted 
to obtain time- resolved spectr oscopy of coron al flares (e.g.. iMacNeice et aZ.I . ll985T ) . Similarly, the 
Hinotori mission (1981 - 1991, iTanakal . Il983h was dedicated mainly to solar flare observations in 
the X -ray band. This was also the scope of the later Yohkoh mission, (1991-2001. lOgawara et al 



199ll ) by means of high resolution X-ray s pectroscopy, adding the monitoring and imaging of the 



hot and flaring corona. Hara et al. ( 19921 ) fo und first indications of plasma at 5-6 MK in active 
regions with the Soft X-ray Telescope (SXT. iTsuneta et~al. , 1991 ). 

Normal incidence optics have been dev eloped in the late 1980s. A n early experiment was the 
Normal Incidence X-ray Telescope (NIXT, iGolub and Herant . 19891 ). which provided a few high 
resolution coronal images in the EUV band. 

More recent spac e missions dedicated t o study the corona have been the Solar and Heliospheric 
Observatory (SoHO. lDomingo et all fl995l) . launched in 1995 a nd still operative, and the Transition 



Region and Coronal Explorer (TRACE, Handy et al. , 19991 ). launched in 1998 and replaced in 



2010 by the Solar Dynamic Observatory instruments. Both SoHO and TRACE were tailored 
to observe the quiet cor ona (below 2 MK). SoHO i mages the whole corona (Extreme ultraviolet 



Imaging Telescope, EIT, iDelaboudiniere et all . Il995l ) and per forms wide band spec troscopy (Solar 



Ultraviolet Measurements of Emitted Radiation, S UMER, Wilhelm et all Il995l ) and (Coronal 
Diagnostic Spectrometer, CDS, Harrison et al. . 19951 ) in the EUV band; TRACE imaged the EUV 
corona with high spatial (0.5 arcsec) and temporal (30 s) resolution. Both SoHO/EIT and TRACE 
are based on normal incidence optics and contain three different EUV filters that provide limited 
thermal diagnostics. 

Thanks to their capabilities, both missions allowed to address finer diagnostics, in particular 
to investigate the fine transverse structuring of coronal loops, both in its geometric and thermal 
components, and the plasma dynamics and the heating mechanisms at a higher level of detail. 
SoHO and TRACE have been complementary in many respects and several studies attempted to 
couple the information from them. 

Among ot her relevant missions, we mention the CORONAS series (Ignatiev et all 19981: Oravevskii 
and Sobelman,l2002J), with instruments like SPectroheliographlc X-Ray Imaging Te lescope (SPIRIT , 



Zhitnik et al. 



_ _ , , Rentgenovsky Spektrometr s Izognutymi Kristalami (ReSIK jSvlwester et all . 

1998) and Solar Photometer in X-rays (SPHINX, ISvlwester et al\ . 120081 ) which have contributed 



to the investigation of coronal loops. 

In late 2006, two other major solar missions started, namel y Hinode (IKosugi et al. 



2003) 



and the Solar TErrestrial Relations Observatory (STEREO) (e.g.. iKaiser et aLl . l2008l ). On-board 
Hinod e, two instruments address particularly the study of coronal loops: the X -Ray Telescope 



(XRT , lGolub et al. . 2007 ) and the Extreme-ultraviolet Imaging Spectrometer (ElS. ICulhane etal 
20071 )). Both these instruments offer considerable improvements on previous missions. The XRT 



has a spatial resolution of about 1 arcsec, a very low scattering and the possibility to switch among 
nine filters and combinations of them. EIS combines well spectral (~ 2 mA), spatial (2") and 
temporal (~ 10 s) resolution to obtain accurate diagnostics of plasma dynamics and density. One 
big achievement of the STEREO mission is that, since it consists of two separate spacecrafts getting 
farther and farther from each other, it allows - through for instance its Sun-Earth Connection 
Coron al and Heliospheric Investigation (SECCHI) p ackage - a first 3D reconstruction of coronal 
loops ( Aschwanden et al. . 20091 : Kramar et al. . 20091 ). 

Hinode and STEREO represent the state-of-art for coronal loop observations, and the coronal 
loop studies including recent findings of these missions will be overviewed. In 2010, the new big 
solar mission Solar Dynamic Observatory (SDO) has been launched and new interesting results 
about the solar corona are expected soon, with instruments like the Atmospheric Imaging Assembly 
(ALA). 
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3 The Observational Framework 



3.1 General properties 

Although coronal loops are often well defined and studied in the EUV band, detected by many 
space mission spectrometers like those on board SoHO and Hinode, and by high resolution imagers 
such as NIXT and TRACE, the bulk of coronal loops is visible in the X-ray band (Figure[T]). Also, 
the peak of the coronal emission measure of active regions - where t he loops are bright est - is 
above 2 MK, which is best observed in X-rays (e.g., Peres et al . 200d : iReale et all l2009af ) . 



Coronal loops are characterized by an arch-like shape that recalls typical magnetic held topol- 
ogy. This shape is replicated over a wide range of dimensions. Referring, for the moment, to the 
soft X-ray band, the main properties of coronal loops are listed in Table [TJ The length of coronal 
loops spans over at least 4 orders of magnitude: bright points (~ 10 8 cm), small Active Region 
(AR) loops (~ 10 9 cm), AR loops (~ 10 10 cm), giant arches (~ 10 11 cm) (Figure [5]). As already 
mentioned, the loops owe their high luminosity and variety to their nature of magnetic flux tubes 
where the plasma is confined and isolated from the surroundings. Magnetized fully-ionized plasma 
conducts thermal energy mostly along the magnetic held lines. Due to the high thermal insulation, 
coronal loops can have different temperatures, from ~ 10 5 K (cool loops), to a few ~ 10 6 (X-ray 
loops), up to a few ~ 10 7 K (flaring loops). A density of the confined plasma below 10 7 — 10 8 cm -3 
can be difficult to detect, while more typical values of bright loops are 10 9 — 10 10 cm -3 in quiescent 
and active regions loops. Flaring loops can be easily a factor 10 denser. The corresponding plasma 
pressure can typically vary between 10~ 3 and 10 dyne cm -2 for non-flaring loops, corresponding 
to confining magnetic fields of the order of 0.1-10 G. One characterizing feature of coronal loops 
is that typically their cross-section is constant along their length above the transition region, at 
variance from the topology of potential magnetic fie lds. There is e vidence that the cross-section 
varies across the transition region, as documented in iGabriell (jl976h . 



Table 1: Typical X-ray coronal loop parameters 


Type 


Length 


Temperature 


Density Pressure 




[10 9 cm] 


[MK] 


[10 9 cm" 3 ] [dyne cm" 2 ] 


Bright points 


0.1-1 


2 


5 3 


Active region 


1-10 


3 


1-10 1-10 


Giant arches 


10-100 


1-2 


0.1-1 0.1 


Flaring loops 


1-10 


> 10 


> 50 > 100 



3.1.1 Classification 



Myriads of loops populate the solar corona and constitute statistical ensembles. Attempts to define 
and classify coronal loops were never easy, and no finally established result exists to-date. Early 
attempts were based on morphologi cal criteria, i.e., bright points, active region loops and large 
scale structures IVaiana et ali (|1973[ ). largely observed with instruments in the X-ray band. In 
addition to such classification, more recently, the observation of loops in different spectral bands 
and the suspicion that the difference lies not only in the band, but also in intrinsic properties, 
have stimulated another classification based on the temperature regime, i.e., cool, warm, hot loops 
(Table [2]). Cool loops are gene rally detected in UV lines at temperatures between 10 5 and 10 6 K. 
They were fir s t add ressed by iFoukall ( 1976 ) and later explored more with SoHO observations 
(jBrekke et al. . 1997 ). Warm loops are well observed by EUV i magers such as So HO/EIT and 
TRACE, and confine plasma at temperature around 1 — 1.5 MK Lenz et al. (|l999l ). Hot loops 
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Figure 1: Images of the same active region, taken in the EUV band with TRACE (top) and in 
the X-ray band with Hinodc/XRT (bottom), on 14 november 2006. The X-ray image shows more 
clearly that the active region is densely populated with coronal loops. 
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are those typically observed in the X-ray band and hot UV lines (e.g. Fexvi), with temperatures 
around or above 2 MK ( Table [J). Th e se are the coronal loops already identified, for instance, in 
the early rocket missions IVaiana et al. (1973). This distinction is not only due to observation with 
different instruments and in different bands, but there are hints that it may be more substantial and 
physical, i.e., there may be two or more classes of loops that may be governed by different regimes 
of physical processes. For instance, the temperature along warm loops appears to be distributed 
uniformly and the density to be higher than that predicted by equilibrium conditions. Does this 
make such loops intrinsically different from hot loops, or is it just the signature that warm loops 
are a transient conditions of hot loops? 



Table 2: Thermal coronal loop classification 



Type 


Temperature [MK] 


Cool 


0.1-1 


Warm 


1-1.5 


Hot 


> 2 







A real progress in the insight into coronal loops is expected from the study of large samples 
of loops or of loop populations. Systematic studies of coronal loops suffer from the problem of 
the sample selection and loop identification, because, for instance, loops in active regions overlap 
along the line sight. Att empts of systematic studies have been performed in th e past on Yohkoh 
and TRACE data (e.g., IPorter and Klimchukl . Il995l : lAschwanden et all l2000l ). A large number 



of loops were analyzed and it was possible to obtain meaningful statistics. However, it is difficult 
to generalize the results because of limited samples and/or selection effects, e.g., best observed 
loops, specific instrument. One basic problem for statistical studies of coronal loops is that it 
is very difficult to define an objective criterion for loop identification. In fact, loops are rarely 
isolated; they coexist with other loops which intersect or even overlap along the line of sight. 
This is especially true in active regions where most of the loops are found. In order to make a 
real progress along this line, we should obtain loop samples and populations selected on totally 
objective and unbiased criteria, which is difficult due to the problems outlined above. Some steps 
are coming in this direction and we will see results in the future. 

3.2 Morphology and fine structuring 
3.2.1 Geometry 

Coronal loops are magnetic structures and might therefore be mapped easily and safely by mapping 
reliably the coronal magnetic field. Unfortunately, it is well-known that it is very difficult to 
measure the magnetic field in the corona, and it can be done only in very special conditions, e.g. 
very strong local field (White et all Il99lh . In some cases it is possible to use coronal seismology 



(first proposed bv lUchida . 1970f) to determ ine t he average magnetic field streng th in an oscillating 



loop first used by iNakariakov et al. ( 1999f ) , and Nakariakov and Ofmanl (|200lh on TRACE loops 



and recently investigated in a number of studies. The accuracy of this method depends on the 
correct detection of the temporally and spatially resolved mode of oscillation, and on the details 
of the loop geometry. 

Since we cannot well determine the coronal magnetic field, coronal loop geometry deserves 
specific analysis. As a good approximation, loops generally have a semicircular shape (Figure [3]). 
The loop aspect, of course, depends on the loop orientation with respect to the line of sight: loops 
with the footpoints on the limb more easily appear as semicircular, as well as loops very inclined 
on the surface near the center of the disk. The assumption of semicircular shape can be useful to 
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Figure 3: Coronal loops have approximately a semicircular shape (image: SDO/AIA, 171 A filter, 
9 July 2010, credit: NASA/SDO). 
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measure the loop length even in the presence of important deformations due to projection effects: 
the de-projected distance of the loop footpoints is the diameter of the arc. However, deviations 
from circularity are rather common and, in general, the detailed analysis of the loop geometry 
is not a trivial task. The accurate determination of the loop geometry is rather important for 
the implications on the magnetic field topology and reconstruction. It is less important for the 
structure and evolution of the confined plasma, which follow the field lines whatever shape they have 
and change little also with moderate changes of the gravity component along th e field lines. First 
works on the accurate determination of the loop geometry date back to the sixties Saito and Billingsl 
(|1964f ) . More specific ones take advantage of stereoscopic views allowed by huge loops during solar 
rotation, with the aid of magnetic field reconstruction methods. These studies find deviations from 
ideal circularity and symmetry, not surprising for such large structures iBerton and Sakurail (|1985l ) . 
The geometry of a specific loop observed with TRACE was measured in the fram ework of a complete 
study including time-dependent hydrodynamic modeling iReale et ali ([2000bl lal) . In that case, the 



discrepancy between the length derived from the distance of the footpoints taken as loop diameter 
and the length measured along the loop itself allowed to assess the loop as elongated. Later, a 
reconstruction of loop geometry was applied to TRACE observation of medium-sized oscillating 
loops, to derive th e prop erties of the oscillations. In this case, a semicircular pattern was applied 



Aschwanden et ali (|2002l ). The importance of the deviations fr om circularity on constraining loop 
oscillations was remarked later iDvmova and Ruderman (j2006h . 

The STEREO mission is actually contributing much to the analysis of loop morphology and 
geometr y, thanks to its u nique capability to observe the Sun simultaneously from different po- 
sitions. iFeng et presented a first stereoscopic reconstruction of the three-dimensional 



shape of magnetic loops in an active region from two different vantage points based on simultane- 
ously recorded STEREO/SECCHI images. They derived parameters of five relatively long loops 
and constraints on the local magnetic field, and found reconstructed loops to be non-planar and 
more curved than field lines extrapolated from SoHO/MDI measurements, probably due to the 
inadequacy of the linear force-free field model used for the extrapolation. A misalignment of 20- 
40 deg betwee n theoretical model and observed loop s has b een quantified from STER EO results 
and discussed ( Sandman et al. , 120091 ; |Pe Rosa et al\ . l2009h lAschwanden et all (|2008bh presented 
triangulations and 3D reconstructions of 30 coronal loops, using the Extreme UltraViolet Imager 
(EUVI) telescopes of both STEREO spacecrafts and deriving a series of loop characteristics, such 
as the loop plane inclination angles, and the coplanarity and circularity. They derived the param- 
eters of seven complete and quite inc lined loops, and found de viations from circularity within 30% 
and less significant from coplanarity. Aschwanden et all (|2009l ) applied a reconstruction method to 
an active region observed with STEREO by combining stereoscopic triangulation of 70 loops and 
addressing mainly density and temperature modeling with a filling factor equivalent to tomographic 
volume rendering. 

Another interesting issue regarding coronal loop geometry is the analysis of the loop cross- 
section, which also provides information about the structure of the coronal magnetic field. 
Yohkoh/SXT allowed for systematic and quantitative studies of loop morphology and showed 
that the cross-section of coronal loops is approximately constant along their length and do not 
increase significantly. More in detail, a systematic analysis of a sample of ten loops showed that 
the loops tend to be only slightly (~ 30%) wider at their midpoints than at their footpoints, while 
for a bipolar field configuration we would expect expansions by factors. One possible explanation 
of this effect is the presence of significant twisting of the magnetic field lines, and therefore the 
development of electric currents and strong deviations from a potential field. The effect might be 
seen either as a twisting of a single loop or as a "braiding" of a bundle of unresolved thin loops. 
At the same time it was found that the variation of width along ea ch loop tends to be m odest, 
implying that the cross section has an approximately circular shape ( Klimchuk et all. 11992). Im- 
plications of these results on the theory of coronal heating are discussed in iKlimchukl (|2000l ). but 
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the conclusion is that none of the current models alone is able to explain all observed properties. 




Figure 4: Magnetic field lines extrapolated from optical magnetogram superposed on a TRACE 
image. Credit: NASA/ESA/LMSAL 



Important information about the internal structuring of coronal loops comes from the joint 
analysis of the photospheric and coronal magnetic field (Figure 2]). An analysis of the magnetic 
field at the footpoints of hot and cool loops showed, among other results, that the magnetic filling 
factor is low er in hot lo ops fO. 05-0.3 out of sunspots) than in warm loops fO.2 — 0.6^ ( Katsukawa 
and Tsuneta" T2005l ). lLopez Fuentes et all (12006) investigated the magnetic structure of loops 



observed with TRACE, applying a linear force-free extrapolation model to SoHO/MDI data and 
comparing the resulting configuration of the magnetic field directly with TRACE images of the 
same region. They confirmed that, whereas the model predicts a significant expansion of the 
magnetic field structures, at least a factor two, from the footpoints to the corona, and a significant 
asymmetry of the structures, because the magnetic field lines starting from the same footpoint 
can diverge to different other end footpoints, these features are not observed: also TRACE loops 
are quite symmetric and their cross-sectio n is constant to a good degree of approximation, as it 
had been already found for Yohkoh loops dKlimchuk et «L,ll992h. The results in Lopez Fuentes 
et al. ( 20061 ) suggest that the tangling of the magnetic flux strands driven by the photospheric 
convection might be very strong, and confirm, therefore, that the magnetic field struct ure is far 
more complicated than it can be modelled even with linear force- free extrapolation. ISchrii vcr 



( 20071 ) studied braiding-induced interchange reconnection of the magnetic field and the width of 
solar coronal loops and showed that loop width observations support the hypothesis that granular 
braiding is countered statistically by frequent coronal reconnections, which in turn explain the 
general absen ce of entangled co ronal field structures in high-resolution observations of the quiescent 
solar corona. IPeForestl ({20071 ) addressed the apparent uniform cross section of bright threadlike 
structures in the corona, long apparent scale height, and the inconsistency between loop densities 
derived by spectral and photometric means. They found that, if coronal loops are interpreted as 
a mixture of diffuse background and very dense, unresolved stranded structures, this requires a 
combination of high plasma density within the structures, which greatly increases the emissivity 
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of the structures, and geometric effects that attenuate the apparent brightness of the feature at 
low altitudes. 

New methods to improve the m orphological analysis of loops are being developed. For instance, 
Dudok de Wit and Aucherej (|2007t ) made multispectral analysis of solar EUV images and explored 
the possibility of separating the different solar structures from a linear combination of images. 
They found source images with more contrast than the original ones. 



3.2.2 Fine structuring 



It has been long claimed (e.g., iGomez et all 119931 ) that coronal loops consist of bundles of thin 
strands, to scales below the current instrumental resolution. The task to investigate this substruc- 
turing is not easy because the thickness of the elementary com ponents may be as small as a few 
km, according to some nanoflare models (e.g., Veksteinl . 2009), and th e measured one goe s down 



to the resolution limit of the most powerful imaging instruments (e.g., IGomez et all Il993f ). First 



limited evidence of fine structuring was the low filling factor inferred for loops observed with NIXT 
(|Di Matteo et all Il999h (see Section |3X2|) . The high spatial resolution achieved by the TRACE 
normal incidence telescope allowed to address the transverse structure of the imaged coronal loops. 
TRACE images visibly show that coronal loops are substructured (Figure [5J . 




Figure 5: TRACE image including a system of coronal loops (9 November 2000, 2 UT). Bundles 
of strands are clearly visible. 



There wer e some early attem pts to study the structure along the single strands in TRACE 
observations ( Testa et all . l2002h . Later, it was shown that in many cases, hot loop structures 



observed with the Yohkoh/SXT are not exactly co-spatial with warm structures observed with the 
SoHO/EIT, which is sensitive in the same bands as TRA CE, nor they cool down to become visible 
to the EIT (jNaeata et all 120031 ; ISchmieder et all |2004 ). 

The detailed morphological comparison of an active region showed that hot loops seen in SXT 
(T > 3 MK) and warm loo ps seen in the SoHO/ EIT 195 A band (T - 1 - 2 MK) are located in 



almost alternating manner ([Nagata et 0/1 12003). The anti-coincidence of the hot and the warm 



loops is conserved for a duration longer than the estimated cooling timescale. These results suggest 
that loops are not isothermal in the transverse direction, rather they have a differential emission 
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measure distribution of modest but finite width that peaks at different temperatures for different 
loops (see Section [33)) . 

In apparent alternative. lWinebarger and Warren (2005) showed that hot monolithic loops visible 



with the Yohkoh/SXT are later resolved as stranded cooler structures with TRACE. However, this 
occurs with a time delay of 1 to 3 hours, much longer than the plasma cooling times, and, therefore, 
correlation can be hardly established between the plasma detected by Yohkoh and that detected 
by TRACE. Direct density measurements of loop plasma made from multi-line observations at 
the solar limb ind i cate n ot very dense plasma and relatively high plasma filling factor (0.2-0.9, 
Ugarte-Urra et all 120051)). This is in th e direction of a moderate structuring of the loops. 



Aschwanden and Nightingale (|2005l ) analyzed specifically and systematically TRACE images 
to search for the thinnest coherent structures that can be resolved with TRACE. They found that 
about 10% of the positions can be fitted with an isothermal model and proposed that, since the 
corresponding structures have a uniform thermal dist ribution, they should be elementary loop 
components, with an average width of about 2000 km. lAschwanden et all (120071) studied statisti- 



cally a large set of coronal loops and found further evidence of elementary loop strands resolved 
by TRACE. 

Other studies based both on models and on analysis of observations independently suggest that 
elementary loop component s should be even finer, with typical cross-sect i ons o f the strands to be 
of the order of 10 - 100 km (jBeveridge et all 120031 ; I Car gill and Klimchukl . 12004 IVeksteinl . l2009f ) . 



3.3 Diagnostics and thermal structuring 

The investigation of the thermal structure of coronal loops is very important for their exhaustive 
physical comprehension and to understand the underlying heating mechanisms. For instance, one 
of the classifications outlined above is based on the loop thermal regime, and, we remark, it is 
debated whether the classification indicates a real physical difference. 

Diagnostics of temperature are not trivial in the corona. No direct measurements are available. 
Since the plasma is optically thin, we receive information integrated on all the plasma column along 
the line of sight. The problem is to separate the distinct contributing thermal components and 
reconstruct the detailed thermal structure along the line of sight. However, even the determination 
of global and average values deserves great attention. 

Moderate diagnostic power is allowed by imaging instruments, by means of multifilter obser- 
vations. Filter ratio m aps provide informat ion about the spatial distribution of temperature and 
emission measure (e.g.. FVaiana et~al. , 19731 ). The emission of an optically thin isothermal plasma 



as measured in a j-th filter passband is: 

I 3 =EM Gj(T) (1) 
where T is the temperature and EM is the emission measure, defined as 

EM = I n 2 dV 
Jv 

where n is the particle density, and V the plasma volume. The ratio R^j of the emission in two 
different filters i,j is then independent of the density, and only a function of the temperature: 

_h _ Gi(T) 

RlJ -J]- G^fj (2) 

The inversion of this relationship provides a value of temperature. 

The limitations of this method are substantial. In particular, one filter ratio value provides one 
temperature value for each pixel; this is a reliable measurement, within experimental errors, as long 
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as the assumption of isothermal plasma approximately holds for the plasma column in the pixel 
along the line of sight. If the plasma is considerably multithermal, the temperature value is an 
average weighted for the instrumental response. Since the response is a highly non-linear function 
of the emitting plasma temperature, it is not trivial to interpret the related maps correctly. In 
addition, it is fundamental to know the instrument response with high precision, in order to avoid 
systematic errors, which propagate dangerously when filter ratios are evaluated. In this sense, 
broadband filters provide robust thermal diagnostics, because they are weakly dependent on the 
details of the atomic physics models, e.g., on the presence of unknown or not well-known spectral 
lines, on the choice of element abundance s. Narrowband filters can show n on-unique dependencies 
of filter ratio values on temperature (e.g., Patsourakos and KlimchulJ . 20071 ). due to the presence of 
several important spectral lines in the bands, but a more general problem can be the bias to detect 
narrow ranges of temperatures forced by the specific instrument characteristics. This problem can 
be important especially when the distr ibution of the emiss ion measure along the line of sight is 
not simple and highly non-linear (e.g., iReale et all 12009 a ). The problem of diagnostics of loop 



plasma from filter ratios, and, more in general, the whole analysis of loop observations, are made 
even more difficult by the invariable presence of other structures intersecting along the line of 
sight. A uniform diffuse background emission also affects the temperature diagnostics, by adding 
systematic offsets which alter the filter ratio values. The task of subtracting this "background 
emission" from the measured emission is non-trivial and can affect seriously the results of the whole 
analysis. This problem emerged dramatically when the analysis of the same large loop structure 
observed with Yohkoh/SXT on the solar li mb led to th r ee different results depending mostly o n 
the different ways to treat the background ( Priest et al . 2000l : Aschwanden . 2001 : lRealel . l2002bh . 
The amount of background depends on the instrument characteristics, such as the passband and 
the point response function: it is most of the si gnal in TRACE UV filterbands, for instance, and its 
subtraction becomes a ve r y delicate issue (e.g., 



Del Zanna and Mason. 2003: Realc and Ciaravella . 



120061 : 1 Aschwanden et al. . 2008a : Terzo and Reale . 2010l ). The problem can be mitigated if one 
analyzes loops as far as possible isolated from other loops, but this is not easy, for instance, in 
active regions. If this is not the case, broadband filters may also include contamination from many 
structures at relatively different temperature and make the analysis of single loops harder. The 
problem of background subtraction in loop analysis has been addressed by several authors, who 
apply different subtraction ranging from simple offset, to emission in nearby pixels or subregions, to 
values interpola ted between the loop sides, to whole images at times when the loop is no longer (or 
not yet) visible dTesta et all 120021: iDel Zanna and Masonl. 120031: TSchmelz et ail 120031: Aschwanden 
and N ightingale, 120051 ; IReale and Ciaravellal . 120061 : lAschwanden et all l2008al ; iTerzo and Reald . 
20101) . 



More accurate diagnostics, although with less time and space resolution, is in principle pro- 
vided by spectrometers and observations in temperature-sensitive spectral lines, which are being 
constantly improved to provide better and better spatial information. Early results from UV spec- 
troscopy alread y recognized the link between transition re gion and coronal loops, for instance, from 



Sky lab mission (jFeldman et a/.l . ll979l ; lMariska et all\198w . Together with background subtraction, 



one major difficulty met by spectroscopic analysis is that, in the UV band, the density of lines is 
so high that they are often blended and, therefore, it is hard to separate the contribution of the 
single lines, especially the weak ones. Fine diagnostics, such as Dopplcr shifts and line broadening, 
can become very tricky in these conditions and results are subject to continuous revisitation and 
warnings from the specialized community. The problem of background subtraction is serious also 
for spectral data, because their lower spatial and temporal resolution determines the presence of 
more structures, and therefore more thermal components, along the line of sight in the same spatial 
element. 

Care should be paid also when assembling information from many spectral lines into a recon- 
struction of the global thermal structure along the line of sight. Methods are well-established 
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Figure 6: Temperature map of an active region obtained from the ratio of two images in different 
broadband filters with Hinode/XRT (12 November 2006, 12 UT). 



(e.g., iGabriel and Jordan! . I1975T) and several ap proaches are available. The so-called method of 



the emission measure loci ( Jordan et all . Il987l ) is able to tell whether plasma is isothermal of 
multithermal along the line of sight (Figure [7]), but less able to add details. Detailed emission 
measure distri butions can b e obtained from differentia l emis sion measure (DEM) reconstruction 
methods (e.g., Brosius et ai , 19961 : iKashvap and Drake) . 1998), but this is an ill-posed mathemat- 
ical problem, and, therefore, results are not unique and are subject to systematic and unknown 
errors. Forward modeling and simulations can be ways to escape from these problems, but they 
require non-trivial computational efforts and programming, and it is not always possible to provide 
accurate confidence levels. All these approaches are constantly improved and probably the best 
way to proceed is to combine different approaches and multiband observations and to finally obtain 
a global consistency. 

In addition to the problems intrinsic to diagnostic techniques, we have to consider that loops 
appear to have different properties in different bands, as mentioned in Section 13.1.11 It is still 
debated whether such differences originate from an observational bias due to the instruments or 
from intrinsic physical differences, or both. In view of this uncertainty, in the following we will 
make a distinction between hot and warm loops, which will generally correspond to loops observed 
in the X-ray (and hot UV lines, e.g., Yohkoh/SXT, Hinode/XRT, SoHO/CDS Fexvi line) and in 
the UV band (e.g., SoHO/EIT, TRACE), respectively. Cool loops are also observed in the UV 
band. The boundary between hot and warm loops is, of course, not sharp, and it is not even clear 
whether they are aspects of the same basic structure, or they really are physically different and 
are heated differently (see also Section l5.3.3p . We will devote attention to the comparison between 
hot and warm loops. 
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3.3.1 Hot loops 

After the pioneering analyses driven by the Sky lab X-ray instruments (Section [5]), Yohkoh/SXT 
allowed to conduct large scale studies on the thermal and structure diagnostics of hot loops, and 
the comparison with other instruments, for instance on-board SoHO, allowed to obtain important 
cross-checks and additional information. Filter ratio maps of flaring loops were shown early after 



the mission launch (jTsuneta et all Il992l ) . 



Systematic measurements of temperature, pressure and len gth of tens of quiescent and active 
region coronal loops were conducted on Yohkoh observations ( Porter and Klimchuk , 19951 ) using 
the filter ratio method. For this sample of loops, selected to be steady and isolated, the lengths 
were measured with assumptions on the loop geometry and ranged in a decade between 5 x 10 9 < 
2L < 5 x 10 10 cm. The temperature measurements were averaged over about half of the loops and 
also ranged in a decade (2 < T < 30 MK), with a mean of about 6 MK. Therefore, it appears as 
a sample of particularly hot loops, although the uncertainties in the hot tail of the distribution 
are very large, probably due to the fiat dependence of the temperature on the filter ratio at high 
temperature. Pressures were derived from the equation of state, after derivation of the density, 
from the emission measure and from the volume inferred from the length and assumptions on 
the loop aspect. They ranged in two decades (0.1 < p < 20 dyne cm~ 2 ). Overall, it was shown 
that the temperature and length of this sample of hot loops are uncorrelated, that pressure varies 
inversely with length (as overall expected for a thermally homogeneous sample from loop scaling 
laws, see Section [4. 1.1 [I . although with a large spread. Such distributions were used as constraints 
on the loop heating through the derivatio n of the dependence of the magnetic field intensity on the 



loop length (Klimchu 



c and Ported . 119951 ) . They also led to accurate analysis of data uncertainties 



( Klimchuk and Garvl . 1995h . Another systemat ic analysis was made on a sample of about 30 bright 



steady Yohkoh loops located in active regions ( Kano and Tsuneta . 19951 ). The temperatures were 
measured after averaging several images and taking the value at the loop top. While this analysis 
confirmed some of the correlations found in IPorter and Klimchukl ( 1995 ) , it found a correlation 
between the loop length and the temperature, and showed deviations from RTV scaling laws 
(Section 14. 1 . 1 [) . It cannot be excluded that correlations between parameters depend on the loop 
sample, as a single scaling law links three parameters. Yohkoh/SXT loops hotter than 3 MK were 
found also in ano ther study, the hottest ones w ith shorter lifetimes (less than few hours) , and often 
exhibiting cusps (jYoshida and Tsunetal . ll996l ). 

A big effort has recently been devoted to the possible detection of hot plasma outsid e of evident 
flares . This would be a conclusive evidence of nanoflaring activity in coronal loop (e.g.. iKlimchukl . 
2006, see Section |4T4"|) . Hinode instruments appear to be able to provide new interesting contribu- 
tions to this topic. The analysis of spectroscopic observations of hot lines in solar active regions 
from Hinodc/EIS allows to construct emission measure distributions in the 1-5 MK temperature 
range, and shows that the distributions are flat or slowly increasing up to approx imately 3 MK 
and then fall off rapidly at higher temperatures (jPatsourakos and KlimchukLT2009f ). Evidence of 
emission from hot lines has been found also in other Hinode/EIS observations, and in particular 
from the analysis of the emission from the Caxvn at 192.858 A, fo rmed near a tem perature of 
6 x 1 6 K, which has be en found in several parts of active regions (|Ko et a&hootf ). Using Fe 



lines, lYoung et all (|2009l ) has shown very accurate density measurements (ss 5%) across an active 



region, with values in the ran ge 8.5 < log (N P J c m 3 ) < 11.0. A smaller density range (from 10 8 5 to 



10 cm 3 has been found bv lWatanabe et all (|2009l ) using the Fe xni line group, although one pair 



has been found to reach the high densit y limi t . Density se nsitive lines have been used to measure 
the filling factor of coronal structures. Dere (2008, 2009) has used spectra and images obtained 
with EIS and comparison with TRACE to determine the volumetric filling factor of bright points. 
The emission measure and bright point widths have been compared with the electron densities and 
with TRACE data. The plasma-filling factor has been found to vary from 3 x 10 -3 to 0.3 with a 



16 



median value of 0.04, which may indicate considerable subresolution structure, or the presence of 
a single completely-filled unresolved loop with subarcsec width. 

Thanks to its multifilter observations, also Hinode/XRT is providing useful information about 
the thermal structure of the bright X-ray corona. Temperature maps derived with combined filter 
ratios show fine structuring to the limit of the instrument resolution and evidence of multithermal 



components (jReale et aZ.1 . 120071 ) . This kind of temperature diagnostics is supported by the evidence 



of warm structures bright in the TRACE images. Observ ations including flare filt ers show evidence 



of a hot component in active regions outside of flares (jSchmelz et all l2009bf ) and data in the 



medium thickness filters appear to constrain better this component of hot plasma a s widespread , 
althou gh minor, and peaking around logT ~ 6.8 — 6.9, with a tail above 10 MK (jReale et all 
20091)). This may support the hypothesis that active regions are heated impulsively. Evidence 



of a persistent although small hot plasma component outside of flare is shown also by RHESSI 



data (jMcTiernanl . l2009h . and the comparison betwee n RHESSI and XRT data seem to support 



this scenario in a consistent way (jReale et all [2009a). The topic of coronal active region heating 
is debated. E videnc e interpreted in the direction of more gradual heating has been obtained by 



Warren et all (|2010f ) (see Section 0~4 



3.3.2 Comparison of hot and warm loops 

Before the SoHO/EIT and TRACE observations, warm loops had been imaged in a similar spectral 
band and with similar optics by the rocket NIXT mission (see Section [2]) . Studies of NIXT loops 
including the comparison with hydrostatic loop models (Section 14.1. lj) poin ted out that brigh t 



spots also visible in Ha band were the footpoints of hot high-pressure loops ( Peres et al . 1994) 



This result was c onfirmed by the comp arison of the temperature structure obtained from Yohkoh 
with NIXT data (|Yoshida et aZj , ll995l i (see also Section [3X3]). 

Another comparison of loops imaged with NIXT and Yohkoh/SXT showed that the compact 
loop structures (length ~ 10 9 cm) have a good general morphologi cal correspondence, whil e larger 



scale NIXT loops (~ 10 10 cm) have no obvious SXT counterpart (|Di Matteo et aZ.I . [l999h . C 



parison with static loop models (see Section I4.1.1j) allowed to derive estimates of the loop filling 
factors, important for the loop fine structure (Section 13. 2. 2j) . In the NIXT band, the filling factor 
of short loops was found to be very low (10~ 3 — 10~ 2 ), but of the order of unity in the SXT band 
and for the largest structure. Information about the loop filling factor was derived also from the 
analysis of simultaneous SoHO and Yohkoh observations of a small solar active region, suggesting 
a volume filling factor decreasing with increas ing density and possi ble differences between emitting 



material in active regions and the quiet Sun ([Griffiths et aLl . l2000l ). 



Some similarity between loops observed in the TRAC E EUV band and Yohkoh X-ray band was 
found out of the core of active region loops iNitta (2000) and interpreted as evidence of loops with 
a broad range of temperatures. Core loops were instead observed only in the X-rays and found to 
be variable, indicating that probably they are not steady. 

The thermal distribution across the loop structures, i.e., along the line of sight, can be investi- 
gated with the analysis of observations in several spectral lines, as obtained, for instance, with the 
SoHO/CDS. Information on the validity of the data analysis and of the loop diagnostics can be 
obtained from the comparison with simultaneous and co-spatial data from imaging instruments. 

Density a nd temperatures in two active regions were accurately determined from SOHO-CDS 
observations ( Mason et al. . 19991 ) and it was confirmed quantitatively that the AR cores are hotter 
than larger loop structures extendi ng above the limb. Fr om the analysis of a single loop observed 
on the solar limb with SoHO/CDS. ISchmelz et al\ (|200ll ) found a bias to obtain flat temperature 
distributions along the loop from ratios of single lines or narrow band filters (TRACE), while a 
careful DEM reconstruction at selected points along the loop is inconsistent with isothermal plasma 
both across and along the loop. A whole line of works started from this study reconsidering and 
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questioning the basic validity of the temperature diagnostics with TRACE and emphasizing the 
importance of the background subtraction, but also the need to obtain ac curate spectral data 
(jSchmelj . |2002| ; iMartens efaH |2002| ; lAschwandenl. |2002|; ISchmelz et a/J.120031). Similar result s but 
different conclusions were reached by lLandi and Landinil (|2004l ) ; lLandi and Feldmanl (|2004l ) who 
analyzed a loop observed with SoHO and, finding it nearly isothermal, considered this evidence as 
real and invoked a non-constant cross-section to explain it. On the other hand, evidence of non- 
uniform temperature along loop s observed with TRACE was also found (|Del Zanna and Masonl . 
20031 : iReale and Ciaravellal . l2006 L emphasizing that the temperature diagnostics with narrow band 
instruments is a delicate issue. 

An interesting debate has focussed on the question whether the loops observed with TRACE 
and CDS have a uniform transverse thermal distribution, i.e., a narrow DEM, or a multi-thermal 
distribution, i.e., a wide DEM which may group together warm and hot loops. Although tackled 
from a different perspective, this question also concerns the fine longitudinal structuring of the 
loop s and of their heating and is th erefore strictly connected to the subject of Sections 13.2.21 and 



4.41 [Del Zanna and Masonl (|2003f ) found a loop detected with TRACE to be isothermal (with 



temperature s below 1 MK) along the line of sight from diagnostics of spectral lines obtained with 
SoHO /CDS. ISchmelz et all (l2005h analyzed another loop on the limb observed with SoHO/CDS, 
with a DEM reconstruction and a careful analysis of background subtraction, and found a multi- 
thermal distribution across th e loop. From the comparison with the isothermal structure of hot 
loops derived from CDS data ( Di Giorgio et aZ.1 . 12003 : Landi and Landini . 2004f ). they concluded 
that there may be two different classes of loops, multi-thermal and isothermal, which they found 
to be confirmed by a systematic inspection of the CDS atlas (Figure [7|). 

Multi band ob servations al low to obtain even more i nformati on and c onstraints. Reale and 
Ciaravella (|2006l ) analyzed a well-defined loop system detected in a time-resolved observation in 
several spectral bands, namely three TRACE UV filters, one Yohkoh/SXT filter, two rasters taken 
with SoHO/CDS in twelve relevant lines (5.4 < logT < 6.4). The data analysis supported a 
coherent scenario across the different bands and instruments, indicating a globally cooling loop 
and the presence of thermal structuring. The study overall indicated that the loop analysis can 
be easily affected by a variety of instrumental biases and uncertainties, for instance due to gross 
background subtraction. The fact that the loop that could be well analyzed across several bands 
and lines is a cooling loop may not be by chance (see end of Section [310]). 

Specific analyses of SoHO spectrometric data have continued to contribute much to the study 
of the loop thermal structure up to recently. A Differential Emission Measure (DEM) analysis of 
coronal loops using a forward- folding technique on SoHO / CDS data has shown different result s 
for two loops, one to be iso thermal and the other to have a broad DEM ( Schmelz et all 2007bl ) . 
Landi and Feldmanl (|2008h have analyzed an extensive active region spectrum observed by the 
SUMER instrument on board SOHO and found that the plasma is made of three distinct isothermal 
components, whose physical properties are similar to coronal hole, quiet-Sun, and active region 
plasmas. 

Hinode/EIS observations of active region loops certify that structures which are clearly dis- 
cerni ble in cooler lines (~ 1 MK) become f uzzy at higher temperatures (~ 2 MK. lTripathi et al 
2009), Figure [8|) , as already pointed out by iBrickhouse and Schmelzl 



Comparative studies of active region loops in the t ransitio n region and the corona Ugarte- 
Urra et al. <|2009h observed with Hinode seem to point out the presence of two dominant loop 
populations, i.e., core multitemperature loops that undergo a continuous process of heating and 
cooling in the full observed temperature range 0.4-2.5 MK shown by the X-Ray Telescope, and 
peripheral loops which evolve mostly in the temperature range 0.4-1.3 MK. 
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Figu re 7: EM-loci p l ots fo r two different loop systems, one showing a multi-thermal structure (top, 
from lSchmelz "aZl . 2005|. reproduced by permission of the AAS), the other an almost isothermal 



one (bottom, from iDi Giorgio et all 120031 ) . The EM is per unit area in the top panel. 
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Figure 8: Loop system observed in several EUV spectral lines with Hinodc/EIS (19 May 2007, 
11:41-16:35 UT). The loops become less and less contrasted, i.e., fuzzier and fuzzier, at higher and 
higher temperature (courtesy D. Tripathi). 



3.3.3 Warm loops 

The TRACE mission opened new intriguing questions because the data showed new features, e.g., 
stranded bright structures mostly localized in active regions, name "the moss" , and because the 
narrow band filters offered some limited thermal diagnostics, but not easy to interpret. Reliable 
temperatures are in fact found in a very narrow range, and many coronal loops are found to be 
isothermal in that range. 

As mentioned in Section 13.3.21 firs t loop diagnostics with normal incidence telescopes were 
obtained from data collected with NIXT [Feres et al\ |l994). The bright spots with Ha counterparts 
were identified with the footpoints of high pressure loops, invisible with NIXT because not sensitive 
to plasma hotter than 1 MK. They have been later addressed as the moss in the TRACE images, 
which undergo th e same effect. Interac tions of moss with underlying chromospheric structures were 
first described bv lBerger et al. | (Il999al ). Comparison of SOHO/CDS and TRACE observations led 



to establish that the plasma responsible for the moss emission has a temperature range of about 
1 MK and is as sociated with hot loops at 1 - 2 MK. with a volume filling factor of order 0.1 ( Fletcher 
and de Pontieu"Il99^~ It was also found that the path along which the emission originates is of 



order 1000 km long. According to an analytical loop model, a filling factor of about 0.1 is in 



agree ment with the hypothesis of moss emission from the legs of 3 MK loops (jMartens et 
2000h . 



As for temperature diagnostics with narrow band filters, loops soon appeared to be mostly 



isothermal with ratios of TRACE filters (jLenz et all . 1 19991 : lAschwanden et all 120001'). Is this a new 



class of loops? Equivalent SoHO /EIT filter ratios provided analogous results (jAschwanden et all , 
1999). This evidence is intriguing and many investigati ons have addressed it (see also Section [3.3.2p . 



From the diagnostic perspective. ISchmelz et a~l\ ( 200lh reconstructed DEM distributions along the 



line of sight from spectral SoHO/CDS data and synthesized EIT count rates from them, which led 
to almost uniform temperatures along the loop, pointing again to an instrumental bias. Weber 
et al. (|2005h confirmed that, provided they are flat, i.e. top-hat-shaped, even broad DEMs along 
the line of sight produce constant TRACE filter ratio values. However, we learn from DEM studies 
made both with spectrometers and from multi- wideband imagers that the DEM of coronal loops is 
most probably neither isothermal nor br oad and flat, instead peaked with com ponents extending 
both to low and high temperatures (e.g.. |Peres et all . I2OO0I : iReale et aZ.I . l2009bl ). The critical point 
becomes the DEM width and it s range of variation. 

Later, ISchmelz et al. I (l2007ah found that even TRACE triple- filter data cannot, in general, con- 
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strain the temperature distribution for plasma in warm loops. On the other hand, Patsourakos and 
Klimchukl2007|) studied the cross-field thermal structure of a sample of coronal loops from triple- 
filter TRACE observations, and found that the observations are compat ible with multi t herma l 
plasma with significant emission measure throughout the range 1-3 MK. ISchmelz et all (|2009ar ) 



used TRACE filter ratios, emission measure loci, and two methods of differential emission measure 
analysis to examine the temperature structure of three different loops. In agreement with previous 
studies, they found both isothermal and multithermal cases. This might not be a contradiction, in 
the view of the presence of at least three possible conditions of warm loops, as discussed at the end 
of this section. iNoglik et all (|2008h compared TRACE to CDS data to measure the temperature 
along a coronal loop in an active region on the solar limb. Their double filter ratio temperature 
analysis technique led to temperatures between 1.0 and 1.3 MK. Emission measure loci from CDS 
lines were consistent with a line-of-sight isothermal structure which increases in temperature from 
~ 1.20 to 1.75 MK along the loop, in contrast with the nearby multithermal background. 

Another puzzling issue, certainly linked to the loop isothermal appearance, is that warm loops 
are of ten diagnos e d to be overdense with respect to the equilibrium values predicted by loop scaling 
laws ( Lenz et al. . 19991 : Winebarger et al. . 2003a . Section [4. 1 . 1)) . To explain both these pieces of 
evidence, several authors claimed that the loops cannot be at equilibrium and that they must 
be filament ed and cooling from a hot ter state, probably continuously subject to heating episodes 
(nanoflares. I Warren et al. . 2002 . 2003 , Sections [4.21 and 14. 4|) . Other authors proposed that p art of 
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the effect might be due to inaccurate background subtraction (|Del Zanna and Mason , 

The Hinode mission is stimulating new analyses of warm coronal structures, mostly based 
on its high quality EIS spectral data. Modeling observations of coronal moss with Hinode/EIS 
confirmed that the moss intensities predicted by steady , uniformly heated loop models are too 



intense relative to the observations (W arren et al. 



2008bl) . A 

and must vary inversely with the loop pressure. Observations of active region loops with EIS 
indicate that isolated coronal loops that are bright in Fe xn generally have very narrow temperature 
distributions (3 x 10 5 K), but are not pro perly isothermal and have a volumetric filling factors of 
approximately 10% (IWarren et al. . 2008a ). 

Schmelz et al. (2008) studied temperatures of loops identified in a TRACE image in three 
density-sensitive line ratios. While emission measure loci plots indicated that the loop plasma 
is not isothermal, a more detailed differential emission measure analysis showed that two broad 
components can reproduce the background-subtracted intensities. They proposed that the two- 
component DEM distribution represents two ensembles of strands, one for each of the loops seen 
in the TRACE image. 

Density diagnostics through density- sensitive line ratio als o led to measure dire ctly density 
values, for instance in active regions fe.g.. lDoschek et al. . 2007bl ). Tripathi et aZ| (l2008h found that 



the hot core of the active region is densest with values as high as 10 cm . The electron density 
estimated in specific regions in the active region moss decreases with increasing temperature. 



The 

density within the moss region was highest at logT = 5.8—6.1, with a value around 10 iu cm~ 3 . 

In a cooler regime (4.15 < logT < 5.45) observed in coordination by SOHO spectrometers and 
imagers, STEREO/EUVI, and Hinode/EIS, active region plasma at the limb has been found to 
cool down from a coronal hole status with temperatures in the 5.6 < logT < 5.9 ranee (X andi 
et aL l2009r i. 

The loop reconstruction analysis described in Aschwanden et al. ( 20091 ) was used mainly for 
density and temperature modeling of the warm loops of an active region observed with STEREO. 
The rendering reduces the problem of background subtraction. Although based on simple model 
assumptions, the derived density and temperature distributions are able to reproduce the total 
observed fluxes within 20%. The modeling extrapolates results quite outside of the range of 
sensitivity of the STEREO EUV filters, anyhow finding emis sion measure distrib utions not very 
different from those obtained from spectroscopic observations ( Brosius et al. , 1996f ) and deviations 
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from hydrostatic values in agreement with other previous studies dLenz et all. 19991: Winebarger 



et aL l2003af ). 

In summary, the current observational framework and loop analysis seems to indicate that for 
a coherent scenario warm loops are manifestations of at least three different loop conditions: i) in 
loops consisting of bundles of thin independently-heated strands, few cooling strands of steady hot 
X-ray loops might be detected as warm overdense loops in the UV band. These warm loops would 
coexist with hot loops and would show a multithermal emission measure distribution fPatsourakos 



and Klimchuk, 200 



3ps and would snow a multithermal emission measure distribution ( ratsourakos 
7l: I Warren et all . l2008al ; iTripathi et all l2009f ) : ii) we might have warm loops as 



an obvious result of a relatively low average heating input in the loop. These loops would be much 
less visible in the X-rays and thus would not be cospatial with hot l oops, and would also be much 
less multithermal (iDi Giorgio et al. , 2003 ; Landi and Landinl . 12004 ; Aschwanden and Nightingale , 



200.4 



Nogli k et al. 



20081) : i ii) wa rm loops might be globally cooling from a status of hot X-ray 
loop (jReale and Ciaravellal . I2006T ) . These loops would also be overdense and cospatial with hot 
loops but with a time shift of the X-ray and UV light curves, i.e., they would be bright in the 
X-rays before they are in the UV band. Also these loops would have a relatively narrow thermal 
distribution along the line of sight. 



3.4 Temporal analysis 

The solar corona is the site of a var iety of trans ient phenomena. Coronal loops sometimes flare 
in active regions (see the review bv iBend . 2008). However, most coronal loops are well-known 
to remain in a steady state for most of their life, much longer than the plasma characteristic 
cooling times (jRosner et al. . Il978l . see Section EXT]). This is taken as an indication that a heating 
mechanism must be on and steady long enough to bring the loop to an equilibrium condition, 
and keep it there for a long time. Nevertheless, the emission of coronal loops is found to vary 
significantly on various timescales, and the temporal analyses of coronal loop data have been used 
to obtain different kinds of information, and as a help to characterize the dynamics and heating 
mechanisms. The time variability of loop emission is generally not trivial to interpret. The problem 
is that the emission is very sensitive to density and less to temperature. Therefore, variations are 
not direct signatures of heating episodes, not even of local compressions, because the plasma is 
free to flow along the magnetic field lines. Variations must therefore be explained in the light of 
the evolution of the whole loop. This typically needs accurate modeling, or, at least, care must be 
paid to many relevant and concurrent effects. 

Another important issue is the band in which we observe. The EUV bands of the normal 
incidence telescopes are quite narrow. Observations are then more sensitive to variations because 
cooling or heating plasma is seen to turn on and off rapidly as it crosses the band sensitivity. 
On the other hand, telescopes in the X-ray band detect hotter plasma which is expected to be 
more sensitive to heating and therefore to vary more promptly, but the bandwidths are large 
and do not take as much advantage of the temperature sensitivity as the narrow bands. Finally, 
spectroscopic observations are, in principle, very sensitive to temperature variations, because they 
observe single lines, but their time cadence is typically low and able to follow variations only on 
large timescales. Time analysis studies can be classified to address two main classes of phenomena: 
temporal variability of steady structures and single transient events, such as flare-like brightenings. 

In spite of limited time coverage, the instrument S-054 on-board Skyl ab already allow ed for early 
studies of variability of hot X-ray loops. Decay times were studied bv iKrieger ( 19781 ) who found 



eviden ce of continued evaporati on of coronal plasma in slowly-decaying structures. iSheelev Jrl 
(|l980h and lHabbal et all (|l985h bund timescales of moderate variability of a few hours over a 
substantial steadiness for observations of active region loops in 2 MK lines such as Fe xv and Si XII. 
Substantial (but non-flaring) temporal variability was reported bv lHaisch etaH (|l988t ) in two active 
region loops observed with SMM in a few relatively hot X-ray lines (~ 5 MK) on time scales of 
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some minutes. Cooler loops (< 1 MK) first studi ed in detail by Foukall ( 1976h (see Section [33 



were found to be more variable and dynamic (e.g., iKopp et all 119851 ) . 



The high time coverage and resolution of Yohkoh triggered studies of brightenings on short 
time scales. The Yohkoh/SXT resolution and dynamic range allowed to study the interaction 
of differently bright hot loops and to show, for instance, that X-ray bright points often involve 
loops considerably large r than the br i ght p oints themselves, and that they vary on timescales 



from minutes to hours ( Strong et al. . 1992). The analysis of a large set (142) of macroscopic 



transient X-ray brightenings indicat ed that they derive from the interaction of multiple loops at 

Some other more specific loop v ariations were also observed, 

19971 ). This was interpreted not 



their footpoints (jShimizu et al 
e.g., the shrinkage of large-scale non- flare loops (jWang et 



199J) 



as an apparent motion, but as a real contraction of coronal loops that brighten due to heating 
at footpoints followed by gradual cooling. Fine-scale motions and brightness variations of the 
emission were found on timescales of 1 minute or less, with dark inclusions corresponding to jets 
of chromospheric plasma seen in the wings of Ha. Such small scale variations are associated with 
the fine structure and dynamic s of the conductively heated upper transition region between the 



solar chromosphere and corona (jBerger et all . [1999b) 



Loop variability was specifically s t udied in several UV spectral lines observed with SoHO/CDS 
for about 3 hours bv lDi Giorgio et al. ( 2003f ). In the hottest lines, within the limited time resolution 
of about 10 min, a few brightenings of a hot loop (~ 2 MK) were detected but they are minor 
perturbations over a steadily high emission level. The observation of the whole life of a cool loop 
(logT ~ 5.3) on a time lapse of a few hours confirmed the highly transient nature of cool loops, 
probably linked to the presence of substantial flows (Section l3.5|) . 

Variability analyses have been conducted also on warm loops present in T RACE data. The 
brightening of a single coronal loop was analyzed in detail bv lReale et all (2000b) in an observation 
of more than 2 hours with a cadence of about 30 s. The loop evolves coherently in the rise phase 
and brightens from the footpoints to the top, allo wing for detailed hvdrodv namic modeling ( Reale 



et al.. l2000al ) (see also Section |4~4| . Active region transient events i.e., short-lived brightenings in 
small-scale lo ops, were obser ved with TRAC E, with a h igh cad ence of 35 s ove r half a n hour ( Seaton 



et al, 200 l[ ). Several brightenings detected over a neutral line in a region of emerging flux were 
interpreted as reconnection events associated with flux emergence, possible EUV counterparts to 
active region tran sient brightenings. Th e fast evolution probably implies high speed flows and high 
coronal densities. IShimojo et all (|2002l ) noticed apparent shrinking and expansion of brightening 
warm loops and proposed heating and cooling of different concentric strands, leading to coronal 
rain visible in the Ha line. Plasma condensations in hot and warm loops were d etected also in 



the a n alysis of line intensity an d velocity in temporal series data from SOHO/CDS (jO'Shea et al 



2007). lAntiochos et all (|2003f ) found no significant variability of the moss regions observed with 



TRA CE. This has been tak en as part of the evidence toward steady coronal heating in active region 



cores (IWarren et all . l2010l Section 14.4 



The analysis of temporal series from various missions has been used, more recently, to investigate 
the possible presence of continuous impulsive heating by nanoflares. The temporal evolution of 
hot coronal loops was studied in data taken with GOES Solar X-ray Imager (SXI), an instrument 



with moderate spatial resolution and spectral band similar to Yohkoh/SXT ([Lopez Fuentes et al 
20071 1 . The durations and characteristic timescales of the emission rise, steady and decay phases 



were found to be much longer than the cooling time and indicate that the loop-averaged heating 
rate increases slowly, reaches a maintenance level, and then decreases slowly (Figure [H]), not in 
contradiction with the early results of Skylab (Section [2]) . This slow evolution is taken as an 
indication of a single heating mechanism operating for the entire lifetime of the loop. If so, the 
timescale of the loop-averaged heating rate might be roughly proportional to the timescale of the 
observed intensity variations. 

Joint TRACE and SOHO/CDS observations allowed to study temperature as a function of time 
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Figure 9: X-ray light curve observed with the SXI telescope on board GOES. The loop lifetime 
is much longer than the characteristic cooling times (courtesy of J. A. Klimchuk and M.C. Lopez 
Fuentes) . 



in active region loops (jCirtain et aliuOOlv . In many locations along the loops, the emission measure 
loci were found consistent with an isothermal structure, but the results also indicated significant 
changes in the loop temperature (between 1 and 2 MK) over the 6 hrs observing period. This was 
interpreted as one more indication of multistranded loops, substructured below the resolution of 
the imager and of the spectrometer. Further support to fine structuring comes from the analysi s 



of the auto-correlation functions in SXT and TRACE loop observations ( Sakamoto et al , 2008) 



The duration of the intensity fluctuations for the hot SXT loops was found to be relatively short 
because of the significant photon noise, but that for the warm TRACE loops agrees well with the 
characteristic cooling timescale. This may support loops to be continuously heated by impulsive 
nanoflares. The energy of nanoflares is estimated to be 10 25 erg for SXT loops and 10 23 erg for 
TRACE loops. The occurrence rate of nanoflares is about 0.4 and 30 nanoflares s _1 in a typical 
hot SXT loop and a typical warm TRACE loop, respectively. 

A recent study on time series has been performed on data taken with the Hinode mission. 
Hinode's Solar Optical Telescope (SOT) magnetograms and high-cadence EIS spectral data were 
taken to study the relationship betwee n chromospheric, t ransition region, and coronal emission 
and the evolution of the magnetic field (jBrooks et al. , 2008| ). The data have allowed to distinguish 



hot, relatively steadily emitting warm coronal loops from isolated transient brightenings and to 
find that they are both associated with highly dynamic magnetic flux regions. Brightenings have 
been typically found in regions of flux collision and cancellation, while warm loops are generally 
rooted in magnetic field regions that are locally unipolar with unmixed flux. The authors suggest 
that the type of heating (transient vs. steady) is related to the structure of the magnetic field, and 
that the heating in transient events may be fundamentally different from that in warm coronal 
loops. 

3.5 Flows 

Diagnosing the presence of significant flows in coronal loops is not an easy task. Apparently moving 
brightness variations may not be a conclusive evidence of plasma motion, since the same effect may 
be produced by the propagation of thermal fronts or waves. Conclusive evidence of plasma motion 
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comes from measurements of Dopplcr shifts in relevant spectral lines. However, the detection 
of significant Doppler shifts requires several conditions to be fulfilled altogether, e.g., significant 
component along the line of sight, amount of moving plasma larger than amount of static plasma, 
plasma motion comparable to typical line broadening effects. 

In general, we can distinguish two main classes of mass bulk motions inside coronal loops: 
siphon flows, due to a pressure difference between the footpoints, and loop filling or draining, due 
to transient heating and subsequent cooling, respectively. Some other evidence of bulk motions, 
such as systematic redshifts in UV lines, has been difficult to interpret. 

Siphon flows have been mainly inv oked to explai n motions in cool loops. The existence of cold 
loops has been known for a long time (jFoukai fl976l ) (see S ection H and SoHO has collected high- 



quality data showing the presence of dynamic cool loops (jBrekke et all 119971 ). A well-identified 



detec tion was found in S oHO/SUMER data, i.e., a small loop showing a supersonic siphon-like 
flow (jTeriaca et all I2004T ) and in SoHO/CDS data (jDi Giorgio et all 120031 ) . 

Redshifts in transition region UV lines have been extensively observed on the solar disk (e.g., 
Doschek et aLl.ll976l:lGebbie et aLl.ll98ll;|Perel.ll982l;lFeldman et aLUl982l:lKlimchukl.ll987l: Rottman 



et al. . ll990l : lBrekkel . Il993l : IPeterl . ll999T ). Some mechanisms have been proposed to explain these red- 
shifts: downward propagating acoustic waves ( Hansteenl. 1993 ). downdrafts d riven by radiatively- 
cooli ng condensations in the solar transition region dReale et al. 1 1996 1 1997bf). nanoflares (T eriaca 



et al, Il999bl ): a conclusive word is still to be given with the better and better definition of the 



observational framework. 

Blue-shift s in the transitio n region are also studied but not necessarily associated with coronal 
loops (e.g., Dere et al . Il986l ). More localized and transient episodes of high velocity outflows, 



name d explosive e vents , have been observed in the transition lines such as C IV , formed at 100,000 K 
(e.g.. |Pere et aZ.I. Il989l: IChae et all Il998bl: Iwinebarger et all Il999l. l2002allbl). However. Teriaca 
et al. (|2002h found indications that such EUV explosive events are not directly relevant in heating 
the corona, are characteristic of structures not obviously connected with the upper corona, and 
have a chromospheric origin . 

Chae et al. (|1998al ) found Doppler shifts increasing and then decreasing with increasing tem- 
perature and explained them by the dominance of emission from plasma flowing downward from 
the upper hot region to the lowe r cool region along flux tubes with varying cross section (a factor 
about 30). iTeriaca et al. (1999a) confirmed these results with the exception of blueshifts at higher 
temperature. 

EUV spectra of coronal loops above an active region show clear evidence of strong dynamical 
activity. In the O v 629 A line, formed at 240,000 K, line-of-sight velocities gr eater than 50 km s" 1 



have been measured with the shift extending over a large fraction of a loop (jBrekke et al\ , \199!v . 



Active region loop structures appear to be e xtremely time variable and dynamic a t transition 
region temperatures, with large Doppler-shifts (|Brekkel . ll999l ). iDi Giorgio et all (|2003l ) showed the 
direct observation and identification of the birth, evolution and cooling of one of such transient 
cool loops, and measured a blue-shifted up-flow all along the loop, probably a one-direction siphon 
flow. 

Winebarger et al. (|2002cf ) analyzed co-aligned TRACE and the SoHO/SUMER observations 
of warm active region loops. Although these loops appear static in the TRACE images, SUMER 
detects line-of-sight flows along the loops of u p to 40 km s -1 . Apparent motions were also detected 



in TRACE images (jWinebarger et all l200ll ). 



In SoHO/EIT high-cadence 304 A images, iDe Groof et al. (2004) analyzed systematic intensity 
variations along an off-limb half loop structure propagating from the top towards the footpoint 
and reported several arguments supporting that these intensity variations are due to flowing/falling 
plasma blobs and not to slow magneto- acoustic waves fSection I4.4|) . This evidence has been the 
object of modeling studies (see Section [473]) . 

The high spectral resolution of Hinodc/EIS is allowing for very detailed studies of persistent 
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Figure 10: Monochromatic (negative) images and dopplergrams (km s 1 ) of an active region 
(NOAA 10926) observed with Hinode/EIS in Fevni, Fexn, Fexv lines. (Courtesy G. Del Zanna) 



loop plasma motions, which are very important to assess whether the loops are to be treated 
more as static or dynamic struc tures as a who l e. Mo re specifically, nonthermal velocities were 



detected in solar active regions ([Doschek et all l2007al ). The largest widths seem to be located 



more in relatively faint zones, some of which also show Doppler o utflows. Doppler flows in active 
region loops observed by Hinode EIS were explicitly addressed bv lDel Zannal (|2008h . who found a 
multifaceted scenario (Figure [TO]) . Persistent redshifts, stronger in cooler lines (about 5-10 km s _1 
in Fe xn and 20 - 30 km s _1 in Fe Vili) , were observed in most loop structures. Persistent blueshifts, 
stronger in the hotter lines (typically 5-20 km s _1 in Fexn and 10-30 km s _1 in Fexv), were 
present in areas of weak emission, in a sharp boundary between the low- lying "hot" 3 MK loops and 
the higher "warm" 1 MK loops. Strong localized outflows (~ 50 km/s ) in a widespread dow nflow 



region were clearly visible in Doppler-shifts maps obtained with EIS (jDoschek et al. . 20081 ). The 
outflows might be tracers o f long loops and/or open magnetic fields. 

Ofman and Wang ( 20081 ) used the high resolution Hinode SOT observations, and detected cool 
plasma flowing in multi-threaded coronal loops with speeds in the range 74-123 km/s. In addition 
to flows, the loops exhibited transverse oscillations. 

Further analysis of coronal plasma motions near footpoints of acti ye region loops sh owed again a 
strong correlation between Doppler velocity and nonthermal velocity (jHara et all 120081 ). Significant 
deviations from a single Gaussian profile were found in the blue wing of t he line profiles f or the 
upflows. These may suggest that there are unresolved high-speed upflows. iTripathi et al. (2009) 
found that an active region was comprised of redshiftcd emissions (downflows) in the core and 
blueshifted emissions (upflows) at the boundary. All these results have to be matched with the 
recent fin ding of extensive blues hifts correlated with spicules upflows and with coronal emission 
intensity (|DePontieu et all l2009h (Section 03}. 
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4 Loop Physics and Modeling 



4.1 Basics 



The basics of loop plasma physics are well established since the 1970s fe.g.. lPriestlll978D . In typical 
coronal conditions, i.e., ratio of thermal and magnetic pressure j3 <C 1, temperature of a few MK, 
density of 10 8 — 10 10 cm -3 , the plasma confined in coronal loops can be assumed as a compressible 
fluid moving and transporting energy only alo ng the magnetic field lines, i.e., along the loop itself 
(e.g.. iRosner et all fl97i IVeseckv et ctZ.I . fl979h . In this configuration, the magnetic field has only 
the role of confining the plasma. It is also customary to assume constant loop cross-section (see 
Section f3 . 2 . 1 p . In these conditions, neglecting gradients across the direction of the field, effects 
of curvature, and transverse waves, the plasma evolution can be described by means of the one- 
dimensional hydrodynamic equations for a compressible fluid, using only the coordinate along the 
loop (Figure fTTjl . 




Figure 11: The plasma confined in a loop can be described with one-dimensional hydrodynamic 
modeling, with a single coordinate (s) along the loop (image: TRACE, 6 November 1999, 2 UT). 



The time-dependent equations of mass, momentum and energy conservation typically include 
the effects of the gravity component along the loop, the radiative losses from an optically thin 
plasma, the plasma thermal conduction, an external heating input, the plasma compressional 
viscosity: 



with p and e defined by: 



and the conductive flux: 
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where n is the hydrogen number density, s the spatial coordinate along the loop, v the plasma 
velocity, tor the mass of hydrogen atom, ji the effective plasma v iscosity, P(T) the radiative losses 
function per unit emission measure (e.g.. iRavmond et al. . 1976[). ffi the frac tional ionization, i.e., 
n c /nn, F c the conductive flux, k the thermal conductivity ( Spitzer . 1 19621 ) . fce the Boltzmann 
constant, and x the hydrogen ionization potential, 
which describes the heat input in the loop. 



H(s,t) is a function of both space and time 



to investigate the physics of coronal loops and of X-ray flares (e.g., Nagai. 1980; Peres et all 1982; 


Doschek et a/.Ul982; Nagai and Emsliel. Il984; 


Fisher et al. 


. 1985allb!ld:lMacNeicel.ll986l:lGan et al. 


1991: iHansteenL Il993: Betta et all 11997c lAntiochos et al. 


. 1999: Muller et all 120031: Ofman and 


Wang.l2002l: Bradshaw and Masonl. 12003: Bradshaw and Cargilll. 20061). 



The concept of numerical loop modeling is to use simulations, first of all, to get insight into 
the physics of coronal loops, i.e., the reaction of confined plasma to external drivers, to describe 
plasma evolution, and to derive predictions to compare with observations. One major target of 
modeling is of course to discriminate between concurrent hypotheses, for instance, regarding the 
heating mechanisms, and to constrain the related parameters. 

The models require to be provided with initial loop conditions and boundary conditions. In view 
of our ignorance of the specific heating mechanisms (see Section the models require to define 
an input heating function, specifying its time-dependence, for instance it can be steady, slowly 
or impulsively changing, and its position in space. The output typically consists of distributions 
of temperature, density and velocity along the loop evolving with time. From simulation results, 
some modelers derive observables, i.e., the plasma emission, which can be compared directly to 
data collected with the telescopes. The model results are, in this case, to be folded with the 
instrumental response. This forward-modeling allows to obtain constraints on model parameters 
and, therefore, qu antitative informati on about the questions to be solved, e.g., the heating rate 
and location (e.g.. iReale et ali 2000a ). 



Loop codes are typically based on finite difference numerical methods. Although they are one- 
dimensional, and therefore typically less demanding than other multi-dimensional codes that study 
systems with more complex geometry, and although they do not include the explicit description 
of the magnetic field, as full MHD codes, loop codes require some special care. One of the main 
difficulties consists in the appropriate resolution of the steep transition region (1-100 km thick) 
between the chromosphere and the corona, which can easily drift up and down depending on the 
dynamics of the event to be simulated. The temperature gradient there is very large due to the 
local balance between the steep temperature dep endence of the thermal conduction and the peak 



of the radiative losses function ([Serio et ali . Il98lh . The density is steep as well so to maintain the 



pressure balance. The transition region can become very narrow during flares. 

Also a fine temporal resolution is extremely important, because the highly efficient thermal 
conduction in a hot magnetized plasma can lead to a very small time step and make execution times 
not so small even nowadays. Another important issue is the necessary presence of a relatively thick, 
cool and dense region under the transition region, i.e., a chromosphere, otherwise the atmosphere 
gets unstable. The main role of the chromosphere in this context is only that of a mass reservoir, 
and therefore, in several codes, it is chosen to treat it as simply as possible, e.g., an isothermal 
inactive layer which neither emits, nor conduct s heat. In other cases , a more accurate description 
is chosen, e.g., to include a detailed model (e.g., Vernazza et aZ.I . fl981 ) and to m aintain a simplified 



radiative emission and a detailed energy balance with an ad hoc heat input ([Peres et al. . Il982l : 



Reale et all l2000af ) 



In recent years, time-dependent loop modeling has revived in the light of the observations with 
SoHO and TRACE for the investigation of the loop dynamics and heating. The upgrade driven by 
the higher quality of the data has consisted in the introduction of more detailed mechanisms for the 
heating input, for the momentum deposition, or others, e.g., the time-dependent ionization and the 
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saturated thermal conduction ( Bradshaw and Cargilll . 20061 : Reale and Orlando! 2008). Some codes 
have been upgraded to include adaptive mesh refinement for better resolutio n in regions of hig h 



gradients, such as in the transition region, or during impulsive events (e.g., iBetta et ali , Il997f h 
Another form of improvement has been the description of loops as collections of thin strands. 
Each strand is a self-standing, isolated and independent atmosphere, to be tre ated exactly as a 



single loop. This approach has been adopted both to descri be loops as static (Reale and Peresl . 



2000) (Figure [12]) and as impulsively heated by nanofiares ([Warren et ali l2002i ). On the same 



line, collect ions of loop models have been applied to describe entire active regions (" Warren and 



Winebarger3o06). 
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Figure 12: E mission in two TRACE filterbands predicted by a model of loop made by several thin 
strands (from Reale and Peresl . 2000l ). 



One limitation of current ID loop models is that they are unable to treat conveniently the ta- 
pering expected going down from the corona to the chromosphere (or expansion upwards) through 
the transition region. This effect can be neglected in many circumstances, but it is becoming 
increasingly important with the finer and finer level of diagnostics allowed by upcoming observa- 
tional data. For instance, the presence of tapering changes considerably the predicted distribution 
of emission measure in the low temperature region (Section HXT| . 

Possible de viations from p ure ID evolution might be driven by intense oscillations or kinks, 
as described in lOfmanl (|2009l ) . The effect of the three-dimensional loop structure should then be 
taken into account to describe the interaction with excited MHD waves (jMcLaughlin and Ofman , 
20081 : IPascoe et ali [20091 : ISelwa and Ofmanl . [2009h . 

However, the real power of ID loop models, that makes them still on the edge, is that they 
fully exploit the property of the confined plasma to evolve as a fluid and practically independently 
of the magnetic field, and that they can include the coronal part, the transition region, and the 
photospheric footpoint in a single model with thermal conduction. In this framework, we may even 
simulate a multi-thread structure only by collecting many single loop models together, still with 
no need to include the description and interaction with the magnetic field. We should however 
be aware that the magnetic confinement of the loop material is not as strong and the thermal 
conduction is not as anis otropic below the coronal part of the loop as it is in the corona. 

Klimchuk et al. (|2008h illustrate a new efficient model of dynamic coronal loops called "Enthalpy- 
Based Thermal Evolution of Loops" (EBTEL), which accurately describes the evolution of the av- 
erage temperature, pressure, and density along a coronal strand with a "0-D" , very fast approach. 
This model is particularly useful for the description of loops as collections of myriads of strands. 
In more detailed modeling, it has been recently shown that non-local thermal conduction may 
lengthen considerably the conduction cooling times and may enhance the chances of observing hot 
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nanoflare-heated plasma ([West et all 120081 ). 

Alternative approaches to single or multiple loop modeling have been developed more recently, 
thanks also to the increasing availability of hig h performance computing syste ms and resources. 
A gl obal "ab initio" approach was presented by Gudiksen and Nordlund (120051) and by Hansteen 



et al. ( 20071 ) (see also lYokovama and Shibatal . l200ll . for the case of a flare model). They model 
a small part of the solar corona in a computational box using a three-dimensional MHD code 
that span the entire solar atmosphere from the upper convection zone to the lower corona. These 
models include non-grey, non-LTE (Local Thermodynamic Equilibrium) radiative transport in the 
photosphere and chromosphere, optically thin radiative losses, as well as magnetic field-aligned 
heat conduction in the transition region and corona. Although such models still cannot resolve 
well fine structures, such as current sheets and the transition region, they certainly represent the 
first important step toward fully self-consistent modeling of the magnetized corona. 



4.1.1 Monolithic (static) loops: scaling laws 

The Skylab mission remarked, and later missions confirmed (Figure [9]), that many X-ray emitting 
coronal loops persist mostly unchange d for a time considera bly longer than their cooling times by 
radiation and/or thermal conduction ( Rosner et al. . 19781 . and references therein). This means 
that, for most of their lives, they can be well descri bed as systems at equilibrium and has been the 
start i ng point for several early theoretical studies (ILandini an d Monsignori Fossil . Il975l : iGabriell . 
19761 : IJordanl . Il976l : IVeseckv et all Il979l : IJordanl . Il980h . iRosner et all (|1978l ) devised a model of 
coronal loops in hydrostatic equilibrium with several realistic simplifying assumptions: symmetry 
with respect to the apex, constant cross section (see Section 13.2. ip , length much shorter than the 
pressure scale height, heat deposited uniformly along the loop, low thermal flux at the base of 
the transition region, i.e., the lower boundary of the model. In these conditions, the pressure is 
uniform all along the loop, which is then described only by the energy balance between the heat 
input and the two main losses mentioned above. From the integration of the equation of energy 
conservation, one obtains the well-known scaling laws: 



and 



T 0)6 = 1.4(pic 
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where Tq^, Lq and -H_3 are the loop maximum temperature To, length L and heating rate per unit 
volume H, measured in units of 10 6 K (MK), 10 9 cm and 10~ 3 erg cm~ 3 s _1 respectively. These 

scaling laws were found in agreement with Skylab data withi n a factor 2. 

Analogous models were developed in the same frame work (ILandini and Monsignori Fossi . 1975 ) 
an d equivalent scaling law s were found independently by Craig etau (1978) and more g eneral ones 
by Hood and Priest ( 19791 ) . They have been derived with a more general formalism by IBrav et HT 



( 199lh . Although scaling laws could explain several observed properties, some features such as 



the emission measure in UV lines and the cool loops above sunspots could not be reproduced 



and, although the laws have been questioned a number of times (e.g., iKano and Tsunetal . Il995l ) 
in front of the acquisition of new data, such as those by Yohkoh and TRACE, they anyhow 
provide a basic physical reference frame to interpret any loop feature. For instance, they provide 
reference equilibrium values even for studies of transient coronal events, the y have allowed to 
constrain that many loop structures observed with TRACE are overdense (e.g.. ILenz et al 1. ll999l : 



Winebarger et al 



status ([Winebarger and Warren 



2003al . Sectior l4".1.2l) and, as such, these loops must be cooling from hotter 
(see Section 13.3.31 and so on. They also are useful for 



density estimates when closed with the equation of state, and for coronal energy budget when 
integrated on relevant volumes and times. 
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Scaling laws have been extended to loops higher than the pressure scale height (ISerio et a 



198ll ) and li mited by the finding that ve ry long loops become unstable ( Wragg and Priestl . Il98lh . 
According to lAntiochos and Nocil (|1986l ). the cool loops belong to a different family and are low- 
lying, and may eventually explain an evidence of excess of emission measure at low temperature. 

The numerical solution of the complete set of hydrostatic equations allowed to obtain detailed 
profiles of the physical quantities along the loop, including the steep transition region. Figure 1131 
shows two examples of solution for different values of heating uniformly distributed along the loop. 
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Figure 13: Distri butions of tempera ture, density and pressure along a hydrostatic loop computed 
from the model o f S erio et ali (|198ll ) for a high pressure loop (AR) and a low pressure one (Empty) 
with heating uniformly distributed along the loop. 



Reale ( 19991 ) and lAschwanden and Nitta (2000) investigated in detail the effect of hydrostatic 
weighting on the loop visibility and on the vertical temperature structure of the solar corona. 
From the comparison of SOHO-CDS observations of active region loops with a static, isobaric loop 
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model (Land ini and L andi. l2002f ). Brko vic et ~d\ <|2002h showed that a classical model is not able 
to reproduce the observations, but ad hoc assumptions might be needed. Further improvements 
of this app roach including flows did not improve the agre ement between the model and the ob- 



servations (jLandi and Landinil . 120041 ) . Using static models. lLandi and Feldmanl (|2004f ) found that 



the loop models overestimate the footpoint emission by orders of magnitude and proposed that 
non-uniformity in the loop cross section, more specifically a significant decrease of the cross section 
near the footpoints, is th e most likely solution to the discrepancy (Section 14. ip . On the same line, 
Winebarger et ali (|2008T ) modelled X-ray loops and EUV moss in an active region core with steady 
uniform heating and found that a filling factor of 8% and loops that expand with height provide the 
best agreement wit h the intensity in t wo X- ray filters, though maintaining still some discrepancies 



with observations. iGontikakis et all (|2008l ) studied the distribution of coronal heating in a solar 



active region using a simple electrodynamic model and attributed the observed small coronal-loop 
width expansion to both the preferential heating of coronal loops of small cross-section variation, 
and the cross-section confinement due to the random electric currents flowing along the loops. 

The strength of scaling laws is certainly their simplicity and their easy and general application, 
even in the wider realm of stellar coronae. However, increasing evidence of dynamically heated, 
fine structured loops is indicating the need for improvements. 



4.1.2 Structured (dynamic) loops 

In the scenario of loops consisting of bundles of thin strands, each strand behaves as an independent 
atmosphere and can be described as an isolated loop itself. If the strands are numerous and heated 
independently, a loop can be globally maintained steady with a sequence of short heat pulses, each 
igniting a single or a few strands (nanoflares) . In this case, although the loop remains steady on 
average for a long time, each strand has a continuously dynamic evolution. The evolution of a loop 
structure under the e ffect of an impulsive heating is well-known an d stud i ed from observations and 
from modeling (e.g. 



rtect or an impulsive neatmg is weii-Known an a stua i ea trom ob servations and 
Nagail. Il980l; IPeres et ali Il982l; ICheng et ali 1 19831 ; iNagai and Emslid . Il984 



Fisher et ali ll985al Jblld: iMacNeicd . Il986: iBetta et aLl^OOll ). since it resembles the evolution o:' 



single coronal flaring loops. We have to mention here that there have been at t empts to model even - 



flaring loops as con sisting of several flaring strands ( Hori et ai , 1997 . 19981 ; Reeves and Warren 



20021 ; IWarrenl . l2006f ). 



The evolution of single coronal loops or single loop strands subject to i mpuls i ve he ating has 
been recently summarized in the context of the diagnostics of stellar flares (|Reald . l2007t ). A heat 
pulse injected in an inactive tenuous strand makes chromospheric plasma expand in the coronal 
section of the strand, and become hot and dense, X-ray bright, coronal plasma. After the end of 
the heat pulse, the plasma begins to cool slowly. In general, the plasma cooling is governed by the 
thermal conduction to the cool chromosphere and by radiation from optically thin conditions. In 
th e following w e outline the evolution of the confined heated plasma into four phases, according 
to iReald (|2007l ) . Figure [14] tracks this evo lution which maps on the path drawn in the density- 
temperature diagram of Figure [15] (see also Ijakimiec et ali Il992l ) . 



Phase I : from the start of the heat pulse to the temperature peak (heating). If the heat pulse 
is triggered in the coronal part of the loop, the heat is efficiently conducted down to the 
much cooler and denser chromosphere. The temperature rapidly increases in the whole loop, 
with a time scale given by the conduction time in a low density plasma (see below). This 
evolution changes only slightly if the heat pulse is deposited near the loop footpoints: the 
conduction front then propagates mainly upwards and on timescales not ve ry different from 
the evaporation time scale s, also because the heat conduction saturates (e.g.. lKlimchukl . 20061 : 



Reale and Orlando . 20081 ). In this case the distinction from Phase II is not clearly marked. 
Phase II : from the temperature peak to the end of the heat pulse (evaporation). The tern 
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Figure 14: Scheme of the evolution of temperature (T, thick solid line), X-ray emission, i.e. the 
light curve (LC, thinner solid line) and density (n, dashed line) in a loop strand ignited by a heat 
pulse. The strand evolution is divided into four phases (I, II, III, IV, see text for further details). 
FromlRealeJ (|2007f l. 



perature settles to the maximum value (To). The chromospheric plasma is strongly heated, 
expands upwards, and fills the loop with much denser plasma. This occurs both if the heating 
is conducted from the highest parts of the corona and if it released directly near the loop foot- 
points. The evaporation is explosive at first, with a timescale given by the isothermal sound 
crossing time (s), since the temperature is approximately uniform in the highly conductive 
corona: 



where m is the average particle mass. After the evaporation front has reached the loop apex, 
the loop continues to fill more gently. The time scale during this more gradual evaporation 
is dictated by the time taken by the cooling rate to balance the heat input rate. 

Phase III : from the end of the heat pulse to the density peak (conductive cooling). When the 
heat pulse stops, the plasma imme diately starts to cool due to the efficient thermal conduction 
(e.g-. TCargill and Klimchukl . Eool , with a time scale (s): 



3n c k B T L 2 = 10.5n c k B L 2 _ ^ m n ^ L l 
2/7kT 7/2 kT 5/2 " T 6 5/2 



1500^ (11) 



where n c (n Cj io) is the particle density (10 10 cm~ 3 ) at the end of the heat pulse, the thermal 
conductivity is n = 9 x 10 (c.g.s. units). Since the plasma is dense, we expect no saturation 
effects in this phase. 

The heat stop time can be generally traced as the time at which the temperature begins to de- 
crease significantly and monotonically. While the conduction cooling dominates, the plasma 
evaporation is still going on and the density increasing. The efficiency of radiation cooling 
increases as well, while the efficiency of conduction cooling decreases with the temperature. 

Phase IV : from the density peak afterwards (Radiative cooling). As soon as the radi ation cooling 



time becomes equal to the conduction cooling time (jCargill and Klimchukl . 120041 ) , the density 
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reaches its maximum, and the loop depletion starts, slowly at first and then progressively 
faster. The pressure begins to decrease inside the loop, and is no longer able to sustain the 
plasma. The radiation becomes the dominant cooling mechanism, with the following time 
scale (s): 

3k B T M 3k B T M _ 3000 ^£ (12) 



n M P{T) bT^n M n M ,9 

where Tm (Tm&) is the temperature at the time of the density maximum (10 7 K), um (^m,9) 
the maximum density (10 9 cm -3 ), P(T) the plasma emissivity per unit emission measure, 
expressed as: 

P(T) = bT a 

with b = 1.5 x 10~ 19 and a = —1/2. The density and the temperature both decrease 
monotonically. 

The presence of significant residual heating could make the decay slower. In single loops, this 
can be di agnosed from the analysis of the slope of the decay path in the density-temperature 
diagram ()Svlwester et all Il993 ; Reale et al. . 1997a ). The free decay has a slope between 



1.5 and 2 in a log density vs log temperature diagram; heated decay path is flatter down to 
a slope ~ 0.5. In non-flaring loops, the effect of residual heating can be mimicked by the 
effect of a strong gravity component, as in long loops perpendicular to the solar surface. The 
depe ndence of the decay slope on the pressure scale height has been first st udied in . Real e 
et al. (|l993l ). and, more recently, in terms of enthalpy flux bv Bradshaw and CargnT ( 2010h . 



As clear from Fig. [15] the path in this phase is totally below, or at most approaches, the 
QSS curve. This means that for a given temperature value the plasma density is higher than 
that expected for an equilibrium loop at that temperature, i.e., the plasma is "overdense". 
Evidence of such overdensity (Section I3.3.3|) has been taken as an important indication of 
steadily pulse-heated loops. 




Figure 15: Scheme of the evolution of pulse-heated loop plasma of Figure [14] in a density- 
temperature diagram (solid line). The four phases are labeled. The locus of the equilibrium loops 
is shown (dashed- dotted line, marked with QSS), as well as the evolution path with an extremely 
long heat pu lse (dashed line) and the corresponding decay path (marked with EQ). Adapted from 
Reald (120071 ). 

This is the evolution of a loop strand ignited by a transient heat pulse. Important properties 
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of the heated plasma can be obtained from the analysis of the evolution after the heating stops, 
i.e., when the plasma cools down. 

Serio et al. (|l99l[ ) derived a global thermodynamic time scale for the pure cooling of heated 
plasma confined in single coronal loops, which has been later refined to be (|Reald . l2007h (s): 



4.8 x 1CT 4 - 



1% 



500- 



Lo 



(13) 



This decay time was obtained assuming that the decay starts from equilibrium conditions, 
i.e., departing from the th e locus of the equilibrium loops with a given length (hereafter QSS 
line, Jakimiec et al. . 1992) in Figure IT5l It is therefore valid as long as there is no considerable 
contribution from the plasma draining to the eenrgy balance. The link between the assumption of 
equilibrium and the plasma evolution is shown in Figure [15] if the heat pulse lasts long enough, 
Phase II extends to the right, and the heated loop asymptotically reaches equilibrium conditions, 
i.e., the horizontal line approaches the QSS line. If the decay starts from equilibrium conditions, 
Phase III is no longer present, and Phase II links directly to Phase IV. Therefore, there is no 
delay between the beginning of the temperature decay and the beginning of the density decay: the 
temperature and the density start to decrease sim ultaneously. Also, the decay will be dominated 



by radiative cooling, except at the very beginning (jSerio et all Il99l[ ). 



The presence of Phase III implies a delay between the temperature peak a nd the density peak 



This delay is often observed both in solar flares (e.g., Sylwcstc r et al 



.e.g., van den Oord et al . 19881 : van den Oord and Mewel . 1989l " Favata et al 

20nd : IStelzer et all l2002h . The presence of this delay, whenever observed, is a signature of a 



1993 ) and i n stellar flares 
20001 : iMaggio et al 



relatively short h eat pu l se, or , in other words, of a decay starting from non-equilibrium conditions. 

According to iReald (|2007l ) , the time taken by the loop to reach equilibrium conditions under 
the action of a constant heating is much longer than the sound crossing time (Equation (|10p). 
which rules the very initial plasma evaporation. As already mentioned, in the late rise phase the 
dynamics become much less important and the interplay between cooling and heating processes 
becomes dominant. The relevant time scale is therefore that reported in Equation (|13p . 

Hydrodynamic simulations confirm that the time required to reach ful l equilibrium scales a s 
the loop cooling time (r s ), and, as shown for instance in Figure [16] (see also I Jakimiec et al. , 1992), 
the time to reach flare steady-state equilibrium is: 



te 



2.3r s 



For t > t eqi the density asymptotically approaches the equilibrium value: 



n = 



T 2 



2a 3 k R L 



= 1.3 x 10 6 



T 2 
L 



where a = 1.4 x 10 3 (c.g.s. units), or 



TlQ = 1.3 



T 2 



(14) 



(15) 



(16) 



If the heat pulse stops before the loop reaches equilibrium conditions, the loop plasma maximum 
densi ty is lower than the value at equilibrium, i.e. the plasma is underdense (jCargill and Klimchukl . 
2004, Section F3~4l . Figure [TBI shows that, after the initial impulsive evaporation on a time scale 
given by Equation (|10[) . the later progressive pressure growth can be approximated with a linear 
trend. Since the temperature is almost constant in this phase, we can approximate that the density 
increases linearly for most of the time. We can then estimate the value of the maximum density 
at the loop apex as: 
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Figure 16: Pressure evolution obtained from a hydrodynamic simulation of a loop strand ignited by 
heat pulses of different duration (0.5, 1, 3 times the loop decay time, see text) and with a continuous 
heati ng. M o st of the rise phase can be reasonably described with a linear trend (dashed lines). 
FromlRealeJ (|2007l ). 



nju«no-^ (17) 

where £m is the time at which the density maximum occurs. 

Phase III ranges between the time at which the heat pulse ends and the time of the density 
maximum. The latter is also the time at whic h the decay path crosses the locus of the equilibrium 
loops (QSS curve). According to Reald ( 2007 ). the temperature Tm at which the maximum density 
occurs is: 

T M = 9 x Kr 4 (n M £) 1/2 (18) 
T Mi6 = 0.9(n M , 9 i 9 ) 1/2 ■ 



or 



We can also derive the duration of Phase III, i.e., the time from the end of the heat pulse to 
the density maximum, as 

Ato-Af « T c lnV> (19) 

where 

and r c (Equation (JXTJ) ) is computed for an appropriate value of the density n c . A good consistency 
with numerical simulations is obtained for n c = (2/3)njvf- 
By combining Equation (fT§j) and Equation (|T7|) we obtain: 
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(20) 



which ranges between 0.2 and 0.8 for typical values of -0 (1.2-2). 

These scalings are related to the evolution of a single strand under the effect of a local heat 
pulse. The strands are below the current instrument spatial resolution and, therefore, we have to 
consider that, if this scenario is valid, we see the envelope of a collection of small scale events. 
The characteristics of the single heat pulses become, therefore, even more difficult to diagnose, 
and the question of their frequency, distribution and size remains open. Also from the point of 
view of the modeling, a detailed description of a multistrand loop implies a much more complex 
and demanding effort. A possible ap proach is to lite r ally b uild a collection of 1-D loop models, 



each with an independent evolution ( Guarrasi et al. . 2010h . One common approach so far has 



been to simulate anyhow the evolution of a single strand, and to assume that, in the presence of 
a multitude of such strands, in the steady state we would see at least one strand at any step of 
the strand evolution. In other words, a collection of nanoflare-heate d strands can be described 
as a whole with the time -average of the evolution of a single strand ( Warren et all 2002, 2003; 



Winebarger et al. . 2003al lbl. see also Section |4~2"|) . Another issue to be explored is whether it is 
possible, and to what extent, to describe a collection of independently-evolving strands as a single 
effective evolving loop. For instance, how does the evolution of a single loop where the heating is 
decreasing slowly compare to the evolution of a collection of independently heated strands, with 
a decreasing rate of ignition? To what extent do we expect coherence and how is it connected to 
the degree of global coherence of the loop heating? Is there any kind of transverse coherence or 
ordered ignition of the strands? It is probably reasonable to describe a multi-stranded loop as a 
single "effective" loop if we can assume that the plasma loses memory of its previous history. This 
certainly occurs in late phases of the evolution when the cooling has been going on for a long time. 



4.2 Fine structuring 

The description and role of fine structuring of coronal loops is certainly a challenge for coronal 
physics, also on the side of modeling, essentially because we have few constraints from observations 
(Section 13. 2. 2j) . One of the first times that the internal structuring of coronal loops have been in- 
voked in a modeling context was for the problem of th e interpreta t ion of the uniform filter ratio 
distribution detected with TRACE along warm loops. iLenz et al. (1999) claimed that standard 
hydrostatic loop m odels with uniform heat ing could not explain such indication of uniform temper- 
ature distribution. Reale and Peres! (120001 ) showed that a uniform filter ratio could be produced by 
the superposition of several t hin strands d e scribe d by standard static models of loops at different 
temperatures. In alternative. lAschwandeij (|200l[ ) showed that long loops heated at the footpoints 
result to be mostly isothermal. The problem with this model is that footpoint-heated loops (with 
heating scale height less than 1/3 of the loop half-length) had been shown to be thermally unstable, 
and therefore they cannot be long-lived, as instead observed. A further alter native is to explain 



observations with steady non-static loops, i.e., with significant flows inside ([Winebarger et al 



l200ll j"2002c. see below - ) . Also this hypothesis does not seem to answer the question (" Patsourakos 
et aL 12004 ) . 

A first step to modeling fine-structured loops is to use multistrand static models. Such models 
show some substan tial inconsistencies wit h observations, e.g., in general they predict too large 
loop cross sections (jReale and Peresl . l2000f ). Such st r ands are conceptually different from the thin 
strands predicted in the nanoflare scenario ( Parker . 1988| ). which imply a highly dynamic evolu- 
tion due to pulsed-heating. The nanoflare scenario is approached in multi-thread loop models, 
convolving the independent hydrodynamic evolution of the plasma confined in each pulse-heated 
strand (see Section 14.31) . These are able to match some more features of the evolution of warm 
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loops observed with TRACE (IWarren et all 12001 120031: IWinebarger et all l2003al ibh. By means 
of detailed hydrodynamic loop modeling, IWarren et ali (|2002l ) found that an ensemble of indepen- 
dently heated strands can be significantly brighter than a static uniformly heated loop and would 
have a flat filter ratio temperature when observed with TRACE. As an extension, time-dependent 
hydrodynamic modeling of an evolving active region loop observed with TRACE showed that a 
loop made as a set of small-scale, impulsively heated strands can generally reproduce the spatial 
and temporal pro perties of the observe d loops, such as a delay between t he appearance of the loop 
in diff erent filters ( Warren et al. , l2003f ) . As an evolution of this approach. IWarren and Winebarger 
( 2006f ) modelled an entire active region for comparison with a SoHO/EIT observation. They made 
potential field extrapolations to compute magnetic field lines and populate these field lines with 
solutions to the hydrostatic loop equations assuming steady, uniform heating. As a result, they 
constrained the link between the heating rate and the magnetic field and size of the structures. 
However, they also found significant discrepancies with the observed EIT emission. 

More recently modeling a loop system as a collection of thin unresolved strand with pulsed 
heating has been used to explain why active reg ions look fuzzie r in ha rder energy bands, i.e. X- 
rays, and/or hotter spectral lines, e.g. Fe XVI ( Tripathi et al . 20091 . Section f3 . 3 . 2 1) . The basic 
reason is that in the dynamic evolution of each strand, the plas ma spends a relatively longer time 



and with a high emission measure at temperature about 3 MK ( Guarrasi et al. . 2010f ) 



Although multistrand models appear much more complex than single loop models and need 
further refinements to match all the observational constraints, as mentioned in Section [4.1.2[ they 
certainly represent an important issue for the future of coronal loop comprehension. 



4.3 Flows 



A generalization of static models of loops (Section l4.1.1j) is represented by models of loops with sta- 
tionary flows, driven by a pressure imbalance between the footpoints (siphon flows). The properties 
of siphon flows have been studied by several authors ( Cargill and Priest! Il980l: IPriestl. Il98ll; iNocil. 
Borrini and Noel ll982l;lAntiochod.[i^^;lThomadJl988l;lMontesinos and ThomadJl989l: Noci 
Thomas and Montesinosl. 1990l: ISpadaro et al'. . 1990l: Thomas and Montesinosl . Il991 



1981 



et al.. 



1989 



Peres et al. 



19921 : Mont esinos and Thomasl . 1l993T r Orla ndo et al .1 (|l995al ) developed a complete 



detailed model of loop siphon flows and used it to explore the space of the solutions and to derive 
an extension of RTV scal ing laws to loops containing subsonic flows. 

Orlando et al. (|1995bl ) explored the conditions for the presence of stationary shocks in critical 
and supersonic siphon flows in coronal loops (Figure I17[) , finding that the shock position depends 
on the volumetric heating rate of the loop, and devising related scaling laws. The presence of 
massive flows may alter the line emission with respec t to static plasma, b ecause of the delay o f 



the moving plasma to settle to ionization equilibrium ( Golub et ali Il989h . ISpadaro et al. ( 1990f ) 



modelled that, even including the effect of ionization non-equilibrium, the UV lines are predicted 
to be blue-shifted by loop models. So non-equilibrium emission from flows cannot explain the 
observed dominant redshifts (Section I3.5[) . The effects of no n-equilibrium of ionizatio n in UV line 
emission from shocked siphon flows are further discussed in lOrlando and Peresl (|l999l ). 

Modeling efforts were devoted in the 1990s to ex plain specifically the extensive evidence of 
redshifted UV lines on the solar disk. iHansteenl (|1993l ) used a hydrodynamic loop model including 



the effects of non-equilibrium of ionization to show that the redshifts might be produced by down- 
ward propagating acoustic waves, possibly stimulate d by nanoflares. By means of two-dimensional 
hydrodynamic simulations, iReale et all (|l996l . ll997b1 ) proposed that the UV redshifts might be due 
to downdrafts driven by radiatively-cooling condensations in the solar transition region. In the 
exploration of the parameter space, they found redshifted components at speeds of several km/s 
for ambient pressure values ranging from those typical of quiet Sun to active regions and predicted 
that redshifts may occur more easily in the higher pressure plasma, typical of active regions. 
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Figu re 17: Exa mple of solutions of a siphon flow loop mod el inclu ding a shock. F rom Orlando and 
Peres dl999|)~ 
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Teriaca et al. ( 1999bl ) explored the idea that the occurrence of nanoflares in a magnetic loop 
around the O vi formation temperature could explain the observed redshift of mid-low transition 
region lines as well as the blueshift observed in low coronal lines (T > 6 x 10 5 K). Observations 
were compared to numerical simulations of the response of the solar atmosphere to an energy 
perturbation of 4 x 10 erg, including non-equilibrium of ionization. Performing an integration 
over the entire period of simulations, they found a redshift in C iv, and a blueshift in O VI and 
Ne VIH, of a few km/s, in rea s onabl e agreement with observations. A similar idea was applied 
by IPatsourakos and Klimchuk (|2006h to make predictions about the presence or absence of non- 
thermal broadening in several spectral lines (e.g., Nevni, Mgx, Fexvn) due to nanoflare-driven 
chromosphcric evaporation. Clearly, the occurrence of such effects in the lines depends considerably 
on the choice of the heat pulse parameters. Therefore, more constraints are needed to make the 
whole model more consistent. In other words, modeling should address specific observations to 
provide more conclusive results. 

Theoretical reasons indicate that flows should be invariably present in coronal loop systems, 
although they may not be necessarily important in the global loop momentum and energy budget. 
For instance, it has been show n that the pr esence of at least moderate flows is necessary to explain 
why we actually see the loops ( Lend . [20041 ) . The loop emission and detection is in fact due to the 
emission from heavy ions, like Fe. In hydrostatic equilibrium conditions, gravitational settlement 
should keep the emitting elements low on the solar surface, and we should not be able to see 
but the loop footpoints. Instead, detailed modeling shows that flows of few km/s are enough to 
drag ions high in the corona by Coulomb coupling and to enhance coronal ion abundances by 
orders of magnitudes. Incidentally, the same modeling shows that, for the same mechanisms, no 
chemical fractionation of coronal plasma with respect to photospheric composition as a function 
of the element First Ionization Potential (FIP) should be present in coronal loops. 

Other studies point instead to the relative unimportance of flows in coronal loops. In particular, 
as already mentioned in Section I3.2.2| by means of steady hydrodynamic loop modeling (i.e., 
assuming equilibrium condit i on, an d therefore dropping the time-dependent terms in Equations ([3]) , 
([!]), (|5|)), IPatsourakos et al. (2004) showe d that flows may not be a ble to explain the evidence of 
isothermal loops, as instead proposed by IWinebarger et all (|2002d ). They found that a heating- 
deposited asymmetrically in a loop is able to drive significant flows in the loop and to enhance 
its density to the levels typically diagnosed from TRACE observations, but it also produces an 
inversion of the temperature distribution and, consequently, a highly structured distribution of the 
relevant filter ratio along the loop, which is not observed. 

Plasma cooling is a mechanism that may drive significant downflows in a loop (e.g., Bradshaw 



mg 

and Cargill, 2005; iBradsha w and Cargill, 2010) . As an ext e nsion of studies on modeling catas- 
trophic cooling in loops (jMuller et aZ.I . |2004 ). iMiiller et all (|2005l ) tried to explain the evidence 
of pr opagating intensity variations observed in the Hen 304 A line with SoHO/EIT f De Groof 
et a/.. I2004I . Section l3~5j) . Two possible driving mechanisms had been proposed: slow magnetoa- 
coustic waves or blobs of cool downfalling plasma. A model of cool downfalling blob triggered in 
a thermally-unstable loop heated at the footpoints gave a qualitative agreement with measured 
speeds and predicted a significant braking in the high-pressure transition region, to be checked in 
future high cadence observations in cool lines. 

Plasma waves have been more recently proposed to have an important role in driving flows 
within loops. Acoustic waves excited by heat pulses at the ch romospheric loop fo otpoints and 



damped by thermal conduction in corona are possible candidates (jTarovan et all . 2005). Even more 



attention receiv e d the propagation of Alfven waves in coronal loops. Hydrodynamic loop modeling 
O'Neill and LI (|2005l ) showed that Alfven waves deposit significant momentum in the plasma, 
and that steady state conditions with significant flows and relatively high density can be reached. 
Analogous results were independently obtai ned with a diffe r ent a p proach: consider i ng a w ind-like 
model to describe a long isothermal loop, iGrappin et all ( 20031 ): iGrappin et al. ( 2005 ) showed 
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that the waves can drive pressure variations along the loop which trigger siphon flow s. Alfv en - 



distu rbances have been recently shown to be amplified by the presence of loop flows (jTarovan 
20091 ). 



As listed above, models predict the development of flows inside coronal loops in a wide variety of 
situations, namely evaporation, draining, siphon flows, waves. The challenge will be to distinguish 
clearly among them and to assess the appropriate weight and importance to them both in the 
spatial and temporal distribution. 



4.4 Heating 

The problem of what heats coronal loops is essentially the problem of coronal heating, and is a 
central issue in the whole solar physics. Although the magnetic origin of coronal heating has been 
well-established since the very first X-ray observations of the corona, the detailed mechanism of 
conversion of magnetic energy into thermal energy is still under intense debate, because a series of 
physi cal effec t s con spire to make the mechanism intrinsically elusive. 

Klimchuk ([20061 ) splits the heating problem into six steps: the identification of the source of 
energy, its conversion into heat, the plasma response to the heating, the spectrum of the emitted 
radiation, the final signature in observables. Outside of analytical approaches, the source and 
conve rsion of energy are typical l y stud ied in detail by means of multi-dimensional full MHD models 
(e.g., Gudiksen and Nordlundl . 2005 ). which, however, are still not able to provide exhaustive 
predictions on the plasma response and complete diagnostics on observables. On the other hand, 
the plasma response is the main target of loop hydrodynamic models, which, instead, are not able 
to treat the heating problem in a self-consistent way (Section |4JJ . 

In the investigation of the source of energy, et al\ (| 19801 ) already pointed out that the 

magnetic field plays an active role in heating the coronal loops. They assumed that the field lines 
are wound continuously by the photospheric convective motions and the generated non-potential 
component is dissipated into heating. Several following studies wer e devoted to t he co nnection 



and scaling of the magnetic energy to the coro nal energy content (jGolub et all Il982t ) and to 



the rate of energy release through reconnection (jGaleev et all Il98lh . The photospheric motions 



are therefore the ultimate energy source and stress the field or generate waves depending on 
whether the tim escale of t h e mot ion is long or short compared to the end-to-end Alfven travel 
time. Following Klimchuk ( 20061 ). dissipation of magnetic stresses can be referred to as Direct 
Current (DC) heating, and dissipation of waves as Alternating Current (AC) heating. 

The question of the conversion of the magnetic energy into heat is also challenging, because 
dissipation is predicted to occur on very small scales or large gradients in the corona by classica l 
theory, unless anomalous dissipation coefficients are invoked. As reviewed by iKlimchukl (|2006l ). 
large gradients may be produced in various ways, involving either magnetic field patterns and their 
evolution, magnetic instabilities such as the kink instability, or velocity pattern, such as turbulence. 
For waves, resonance absorption and phase mixing may be additional viable mechanisms. 

The problem of plasma response to heating has been kept historically well separated from the 
primary heating o rigin, although some attempts have been made to couple them. For instance, in 
Reale et all |2005) the time-dependent distribution of energy dissipation along the loop obtained 



from a hybrid shell model was used as heating input of a time-dependent hydrodynamic loop model 
(see below). A similar concept w as applied to search for signatures of turbulent heating in UV 



spectral lines ( Parenti et al . 20061 ) 



As already mentioned, studies using steady-state or time-dependent purely hydrodynamic loop 
modeling have addressed primarily the plasma response to heating, and also its radiative emission 
and the detai led comparis o n with observations. A forward-modeling including all these steps was 
performed bv lReale et ail (|2000bl lah. They first analyzed a TRACE observation of a brightening 
coronal loop (see also Section [3^4]) . The analysis was used to set up the parameters for the forward 
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modeling, and to run loop hydrodynamic simulations with various assumptions on the heating 
location and time dependence. The comparison of the TRACE emission predicted by the simula- 
tions with the measured one constrained the heat pulse to be short, much less than the observed 
loop rise phase, and intense, appropriate for a 3 MK loop, and its location to be probably midway 
between the apex and one of the footpoints. 

The investigation of the heating mechanisms through the plasma response is made difficult by 
a variety of reasons. For instance, the problem of background subtraction can be crucial in the 
comparison with observations, as shown by the three analyses of the same large loop structure 
observed with Yohkoh/SXT on the solar limb, mentioned in Section 13.31 More specifically. Priest 



et al. (|2000h tried to deduce the form of the heating from Yohkoh observed temperature profiles 
and found that a uniform heating best describes the data, if the t emperature is obtai ned from the 
ratio of the total filter intensities, with no background subtraction. lAschwandenl ([2001) splitted the 
measured emission i nto two comp onents and found a better agreement with heating deposited at 
the loop footpoints. iReald (l2002bh revisited the analysis of the same loop system, considering con- 
ventional hydrostatic single-loop models and accounting accurately for an unstructured background 
contribution. With forward- modeling, i.e., synthesizing from the model observable quantities to be 
compared directly with the data, background-subtracted data are fitted with acceptable statistical 
significance by a model of relatively hot loop (~ 3.7 MK) heated at the apex, but it was pointed 
out the importance of background subtraction and the necessity of more specialized observations 
to address this question. More diagnostic techniqu es to compare models with observations were 
proposed afterwards (e.g., Landi and Landini . 20051 ). 

Independently of the adopted numerical or theoretical tool, many studies have been addressing 
the mechanisms of coronal loop heating clearly distinguishing between the two mai n classes, i.e. , 
DC heating through moderate and freq uent explosive events, nam ed nanoflares (e.g. . iParkerl . Il988l ) 
and AC heating via Alfven waves (e.g.. iLitwin and Rosnerl . [l998h . 



4.4.1 DC heating 

Heating by nanoflares has a long story as a possible c andidate to explain the heating of the solar 
corona, and, in particular, of the coronal loops (e.g.. Peres et all 19931: Cargili 19931; Kopp and 



Poletto. Il993l; IShimizul. Il995l: J Judge et all Il998l; iMitra-Kraev and Bend. l200ll: Katsukawa and 
Tsuneta. l200ll; IWarren et all |2002|. 120031; ISpadaro et all 120031; ICargill and Klimchukl . Il997l. 12004; 
Miiller et all 120041 : iTesta et all 120051 : iReale et all 120051 : IVeksteinl . l2009h . ~ 

More specifically, ICarg illl (|l994h provided detailed predictions from a model of loops made of 
thousands of nanoflare-heated strands. In particular, whereas the l oop total emission m easure 
distr ibution should steepen above the canonical T 15 (Ijordanl. 1980l; Orlando et al. k 2000l: Peres 



et aZ.. l200ll ) dependence for temperature above 1 MK. Moreover, it was stressed the importance 
of the dependence of effects such as the plasma dynamics (filling and draining) on the l oop fi lling 
factor driven by the elemental heat pulse size (Section 14.1. 2 1) . In Cargili and Klimchuk ( 1997h the 
nanoflare model was applied to the heating of coronal loops observed by Yohkoh. A good match 
was found only for hot (4 MK) loops, with filling factors less than 0.1, so that it was hypothesized 
the existence of two distinct classes of hot loops. 

Although there is evidence of intermittent heating episodes, it has been questione d whether and 



Aschwandcn 



to wh at extent nanoflares are able to provide enough energy to heat the corona (e.g 
1999). On the other hand, loop models with nanoflares, and, in particular, those considering a 
prescribed random time distribution of the pulses deposited at the footpoints of multi-stranded 
loops have bee n able to explain several f eatures of loop observations, for instance, of warm loops 
from TRACE (|Warren et all |2002| . l2003h . (see Section [3JL2]) . 

Hydrodynamic loop modeling showed also that differen t distributions of the heat pulses along 
the loop have limited effects on the observable quantities ([Patsourakos and Klimchukl . I2005K be- 
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cause most of the differences occur at the beginning of the heat deposition, when the emission 
measure is low, wh ile late r the loop loses memory of the heat distributio n (see also Winebarger 



and Warren" l2004l ). IPatsourakos and Klimchukl (|2008f T applied both static and impulsive models to 
solar active regions and showed that the latter ones are able to simultaneously reproduce EUV and 
SXR loops in active regions, a nd to predict radial inte n sity v ariations consistent with the localized 
core and extended emissions. I Car gill and Klimchukl (j2004h showed with a semi-analytical loop 
model that the cycle of loop hea ting/cooling naturally lea ds to hot-underdense/warm-overdense 



loop (Section [4.1. 2[) , as observed (jWinebarger et all l2003bl . Section [3.3. 3[) , and that the width of 



the DEM of a nanoflare-heated loop can depend on the number of strands which compose the loop: 
a relatively fiat DEM or a peaked (isothermal) DEM are obtained with strands of diameter about 
15 km or about 150 km, respectively. This is of relevance for the diagnostics both of the loop fine 

As a further improvement, 



structure (Section 13.2.21) and of the DEM reconstruction (Section I3.3|) 

Warren and Winebargerl ( 20071 ) added an impulsive heating model to the simulation of an entire 
active region and found that it is possible to reproduce the total observed soft X-ray emission in 
all of the Yohkoh/SXT filters. However, once again, at EUV wavelengths the agreement between 
the simulation and the observation is only partial. 

(2005) 



Testa et al. 



Nanoflares have been studied also in the framework of stellar coronae. 
showed that intermittent heating by relatively intense nanoflares deposited at the loop footpoints 
make the loop stable on long time scales (loops continuously heated at the footpoints are unstable) , 
and, on the other hand, produces a well-defined peak in the average DEM of the loop, si milar to 
that d erived from the DEM reconstruction of active stars, and also to those shown in ICargilll 
(1994). Therefore, this is an alternative way to obtain a steep temperature dependence of the loop 
emission measure distribution in the low temperature range. 

An alternative approach to study nanoflarc heating is to analyze intensity fluctuations (Shimizu 
and Tsuneta, 1997; Vckstcin and Katsukawa, 2000; Katsukawa and Tsuneta, 2001; Vckstein and 



Jain, 120031 ) and to derive their occurrence distribution (jSakamoto et all . l2008l . 120091 ) . From the 
width of the distributions and autocorrelation functions, it has been suggested that nanoflare signa- 
tures are more easily found in observations of warm TRACE loops than of hot Yohkoh/SXT loops. 
It is to be investigated whether the results change after relaxing the assumption of temperature- 
independent distribution widths. Also other variability analysis of TR ACE observations was fo und 
able to put constraints on loop heating. In particular, according to Antiochos et all (|2003l ). in 
TRACE observations, the lack of observable warm loops and of significant variations in the moss 
regions implies that the heating in the hot moss loops should not be truly flare-like, but instead 
quasi-steady and that the heating magni tude is only weakly varying. Further evidence in this 
direction has been found more recently by I Warren et al. ( 2010h . 

An analogous approach is to analyze the intensity distributions. The distribution of impulsive 
events vs their number in the solar and stellar corona is typically described with a power law. The 
slope of the p ower law is a c ritical parameter to establish weather such events are able to heat the 
solar corona ((Hudson, 1991 ). In particular, a power law index of 2 is the critical value above or 
below which flare-l ike events may be able or unable, respectively, to power the whole corona (e.g. 



Aschwanden . 1999f ). Unfortunately, due to the faintness of the events, the distribution of weak 



events is pa rticularly difficult to deriv e and might even be separate from that of proper flares and 
microflares. IParenti and Youngl ( 20081 ) used a hydrodynamic model to simulate the UV emission of 
a loop system heated by nanoflares on small, spatially unresolved scales. The simulations confirm 
previous results that several spectral lin es have an inte nsity distribution that follows a power-law, 
in a similar way to the heating function (|Hudsonl . ll99lh . However, only the high temperature lines 
best preserve the heating function's power law index (especially Fexix). 
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4.4.2 AC heating 

Loop oscillations, modes and wave propagation deserve a review by themselves, and are outside 
of the scope of the present one. Here we account for some aspects which are relevant for the 
loop heat ing. A recent review of coronal waves a nd oscillatio ns can be found in Nakariakov and 
Verwichte"l2005|). New obser vations from SDO AIA provide ample evidence of wave activity in 
the solar corona (Title, 2010t ). These observations are currently the subject of intensive analysis 



and will be report ed on in the futur e. 

As reviewed bv lKlimchukl(|2006h . MHD waves of many types are generated in the photosphere, 
e.g. acoustic, Alfven, fast and slow magnetosonic waves. Propagating upwards, the waves may 
transfer energy to the coronal part of the loops. The question is what fraction of the wave flux 
is able to pass through the very steep density and temperature gradients in the transition region. 
Acoustic and slow-mod e waves form shocks and are stron gly damped, fast-mode waves are strongly 
refracted and ref l ected dNarain and Ulmschneider , 1996f ) . 

Ionson ( 19781 . KL983 . Il983h devis ed an LRC equi valent circuit to show the potential importance 
of AC processes to heat the corona. IHollwesJ (|l984T ) used a dissipation length formalism to propose 
resonance absorption of Alfven waves as a potential coronal heating mechanism. A loop may be 
considered as a high-quality resonance cavity for hydromagnetic waves. Turbulent photospheric 
motions can excite small-scale waves. Most Alfven waves are strongly reflected in the chromosphere 
and transition region, where the Alfven speed increases dramatically with height. Significant 
transmission is possible only within na rrow frequency bands centered on discrete values where 
loop resonance conditions are satisfied ( Hollwegl . Il98li . Il984l : llonsonl . fl982l ). The waves resonate 
as a global mode and dissipate efficiently when their frequency is near the lo cal Alfven waves 
frequency uja ~ 2-kva/ L. By solving the linearized MHD equations iDavilal (|l987f ) showed that this 
mechanism can potentially hea t the corona, as further supp orted by nume rical solution of MHD 
equations for low beta plasma ( Steinolfson and Davilal . 1993 ) , and although IParkerl ( 1991 ) argued 
that Alfven waves are difficult to be generated by solar convection. 

Evidence for photospheric Al fven waves wa s obtained from magnetic and velocity fluctuations in 
regio ns of strong magnetic field (IUlrichl.ll996l) and from granular motions in the quiet Sun (Muller 



et al., 1 1994T ) with fluxes of the order of 10 7 erg cm 2 s 1 , which might contribute to heating if 
transmitted efficiently to the corona. 

Hollweg (|l985h estimated that enough flux may pass through the base of long (> 10 10 cm) ac- 
tive region loops to provide their heating, but shorter loops are a problem, since they have higher 
resonance frequencies and the photospheric powe r spec trum is believed to decrease exponentially 
with frequency in this range. Litwin and Rosner (|l998h suggested that short loops may transmit 
waves with low frequencies, as long as the field is sufficiently twisted. iHollweg and Yangj (|l988h 
proposed that Alfven resonance can pump energy out of the surface wave into thin layers surround- 
ing the resonant field lines and that the energy can be distributed by an eddy viscosity throughout 
large portions of coronal active region loops. 

Waves may be generated directly in the corona, and evidence for their presence was found 
(e.g.. [Nakariakov et all 1 19991; Aschwanden et all 19991; Berghmans and Clettei. 19991: De Moor- 
tel et al.. I2002T). It is unclear whether coronal wav es carry a sufficient energy flux to heat the 



plasma (|Tomczvk et all 120071 ). lOfman et al\ (|l995h studied the dependence on the wavenumber 



for comparison with observations of loop oscillations and found partial agreement with velocity 
amplitudes measured from non-thermal broadening of soft X-ray lines. The observed nonthermal 
broadening of transition region and coronal spectral lines implies fluxes that may be sufficient to 
heat both the quiet Sun and active regions, but it is unclear whether the waves are efficiently 



dissipated (jPorter et al\ . \1994 ). Furthermore, the non thermal line broadening cou l d be produced 



by u nresolved loop flows that are unrelated to waves (Patsourakos and Klimchukl. 200(jf). Ofman 



et al. (|l998l ) included inhomogeneous density structure and found that a broadband wave spectrum 
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becomes necessary for efficient resonance and that it fragments the loop into many density layers 
that resemble multistran d concept. The heat deposition by the resonance of Alfven waves in a 
loop was inves tigated bv lO'Neill and Li ( 2005 ) By assuming a functional form first proposed by 
HollweeJ (1986), hydrodynamic loop modeling showed that, depending on the model parameters, 
heating by Alfven waves leads to different classes of loop solutions, such as the isothermal cool 
loops indicated by TRACE, or the hot loops observed with Yohkoh/SXT. Specific diagnostics are 
still to be defined for the comparison with observations. 

Efficient wave dissipation may be allowed by enh anced dissipation coeffici ents inferred from fast 
damping of flaring loop oscillations in the corona ( Nakariakov et al. , 1999f ). but the same effect 
may also favor efficient magnetic reconnection in nanoflares. Alfven waves required for resonant 
absorption are relatively high frequency waves. Evide nce for lower frequency Alfven waves has 



been found in the chromosphere with the Hinode SOT (|De Pontieu et all 120071 ) . Such waves may 
supply energy in the corona even outside of resonance with different mechanisms to be explored with 
mo deling. Among dissipatio n mech anisms phase mixing with enhance d resistivity was suggeste d 
by lOfman and Aschwanden (|2002T ) and supported by the analysis of lOfman and Wand (|2008h . 
Also multistrand structure has been r ecognized to b e important in possible wave dissipation and 
loop twisting, as recently modelled bv lOfmanl (l2009t ). 

Intensity disturbances propagating along active region loops at spe eds above 100 km / s wer e 
detected with TRACE and interpreted as slow magnetoacoustic waves ( Nakariakov et «l1 .|2000). 
These waves probably originate from the underlying oscillations, i.e., the 3-minute 
chromosphcric/ transition- region oscillations in sunspots and the 5-minute solar global oscillations 
(p- modes). Slow magnetosonic wa ves might be g ood c andidates as coronal heating sources accord- 
ing to quite a detailed model by Belien et al. (1999), including the effect of chromosphere and 
transition region and of the radiative losses in the corona. Such waves might be generated directly 
from upward propagating Alfven waves. Contrary conclusion s, in favor of fas t magn etosonic waves, 



have been also obtained, but with much simpler modeling (jPekunlu et all 120011 ) . Slow magne 



tosonic waves with periods of about 5 minutes have been more recently detected in the transition 
region and coronal emission lines by Hinode/EIS at the foo tpoint of a coronal loop rooted at plage, 
but found to carry not enough energy to heat the corona (I Wang et aZ.1 . 12009^ . 

I nves tigation of AC heatin g has been made also throug h comparison with DC heating. Antolin 
et al. (|2008l ) compared observational signatures of coronal heating by Alfven waves and nanoflares 
using two coronal loop models and found that Hinode XRT intensity histograms display power- 
law distributions whose indices differ considerably, to be checked against observations. Lundquist 
et al. (|2008al lbl) applied a method for predicting active region coronal emissions using magnetic 
field measurements and a chosen heating relationship to 10 active regions. With their forward- 
modeling, they found a volumetric coronal heating rate proportional to magnetic field and inversely 
proportional to field-line loop length, which seems to point to, although not conclu sively, the steady- 
state scaling of two heating mechanisms: van Ballegooijen's current laye rs theory ( van Ballegooijenl . 
198G), taken in the AC limit, and Parker's critical angle mechanism ([Parker . 1 19881 ) . in the case 
where the angle of misalignment is a twist angle. 



4.4.3 Modeling including the magnetic field 

Some studies have investigated loop heating by including also the magnetic field in the analysis and 
modeling, and have tried to discriminate between different mechanisms using global approaches 
and scalings inferred f rom modeling. Based o n a previous study of t he plasma parameters and the 
magnetic flux density ( Mandrini et al . 2000l ). Demoulin et al. I (|2003T) derived the dependence of the 
mean coronal heating rate on the magnetic flux density from the analysis of an active region. By 
using the scaling laws of coronal loops, they found that models based on the dissipation of stressed, 
current-carrying magnetic fields are in better agreement with the observations than models that 
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attribute coronal hea ting to the dissipation of MHD waves injected at the base of the corona. 
Schriiver et all (J200J) considered a similar approach applied to the whole corona, by populating 
magnetic field lines taken from observed magnetograms with quasi-static loop atmospheres and 
obtaining the best match to X-ray and EUV observation with a heating that scales as it would be 
expected from DC reconnection at tangential discontinuities. These approaches will certainly be 
very useful when they will provide more detailed predictions and constraints. 

Recent modeling has been able to explain the ignition of warm loops from primary energy 
release mechanisms, although it is remains unclear how the same mechanisms could produce hot 
loops. A large scale approach (see also Sectio n HZD is by "ab initio" m odeling, i.e., with full MHD 
modeling of an entire coronal region (jGudiksen and Nor dlundl . l2005h . Observed solar granular 
velocity pattern, a potential extrapolation of a SOHO /MDI magnetogram, and a standard stratified 
atmosphere were used as initial conditions. The simulation showed that, at steady state, the 
magnetic field is able to dissipate (3 — 4) x 10 6 erg cm" 2 s" 1 in a highly intermittent corona, at an 
average temperature of ~ 10 6 K, adequate to reproduce typical warm loop populations observed 
in TRACE images. Warm loops were also obtained with time-dependent loop modeling including 
the intermittent magnetic dissipation in MHD turbulence due to loop footpoint motions fReale 



aL~E)05). 



et 



The dissipation rate along a loop predicted with a hybrid-shell model (Ni gro et 



2004) was used as heat ing input ( s ee Eq uation ([5])) in a proper time-dependent loop model, the 



Palermo- Harvard code ([Peres et aUll982l ). It was shown that the most intense nanofiares excited in 
an ambient magnetic field of about 10 G can produce warm loops with temperatures of 1 - 1.5 MK 
in the corona of a 30,000 km long loop. 

More recently. iRappazzo et al. ( 2007 ') used reduced MHD (rMHD) to identify MHD anisotropic 
turbulence as the physical mechanism responsible for the transport of energy from the large scales, 
where energy is injected by photospheric motions, to the small scales, where it is dissipated. Strong 
turbulence was found for weak axial magnetic f ields and long loop s. The predicted heating r ate is 
appropriate for warm loops, in agreement with Reale et a l. (2005). Buchlin and Vellil (2007) also 
used shell models of rMHD turbulence to analyze the case of a coronal loop heated by photospheric 
turbulence and found that the Alfven waves interact non-linearly and form turbulent spectra. They 
derived an intermittent heating function, on avera ge able to sust a in the corona and proportional 
to the aspect ratio of the loop to the ~ 1.5 power. Buchlin et al\ ( 2007 ) added in the modeling a 
profile of density and/or magnetic field along the loop, showing that differences are found in the 
heat deposition, in particular in the low part of the loop. 

(2008) 



There are new efforts to include magnetic effects in the loop modeling. lHavnes et al 



studied observational properties of a kink unstable corona l loop, using a fluid code and finding 
potentially observable density effects. Browning et al. ( 20081 ) studied coronal heating by nanofiares 
triggered by a kink instability using three-dimensional magnetohydrodynamic numerical simula- 
tions of energy release for a cylindrical coronal loop model. Magnetic energy is dissipated, leading 
to large or small heating events according to the initial current profile. 

Interesting perspectives ar e developing from mo dels in which self-or ganized criticality trig- 
gers loop coronal heating. For lUzdenskvl (|2007f ) and ICassak et al. (120081 ) coronal heating is self- 
regulating and keeps the coronal plasma roughly marginally collisionless. In the long run, the 
coronal heating process may be represented by repeating cycles that consist of fast reconnection 
events (i.e., nanofiares), followed by rapid evaporation episodes, followed by relatively long periods 
(~ 1 hr) duri ng which magnetic stresses build up and the plasma simultaneously cools down and 
precipitates. iMorales and Charbonneau ( 2008 ) proposed an avalanche model for solar flares, based 
on an idealized representation of a coronal loop as a bundle of magnetic flux strands wrapping 
around one another. The system is driven by random deformation of the strands, and a form of 
reconnection is assumed to take place when the angle subtended by two strands crossing at the 
same lattice site exceeds some preset threshold. For a generic coronal loop of length 10 10 cm and 
diameter 10 8 cm the mechanism leads to flare energies ranging between 10 23 and 10 29 erg, for an 
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instability threshold angle of 11 degrees between contiguous magnetic flux strands. 



4.4.4 New hints 

Given the difficulty to find a conclusive answer about the heating of coronal loops, even in the 
presence of considerable observational and theoretical efforts, there has been recently the attempt 
to propose different or radically alternative scenarios. For instance, it has been suggested that 
warm and hot lo ops may be heated by different mechanisms, impulsive the former, much more 
steady the latter ( Warren et aZ.L [20 1 dh . New evidence from Hinode satellite indicates the presence 



of significant upflows in the fo r m of widespread spicules which have correspondence in coronal 
observations ( DePontieu et al , 20091 ). This evidence suggests that the interaction between the 



chromosphere and the corona in the heating processes might be imp ortant and may receive support 
also by some modeling approaches (jGudiksen and Nordlund . [2005) . These scenarios are intriguing 
but need additional investigation and support both from the observational and theoretical point 
of view. 
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5 Stellar Coronal Loops 



Non-solar X-ray missions since Einstein and European X-ray Observatory SATellite EXOSAT have 
established that most other stars have a confined corona, often much more active that the solar 
one. The level of activity is ruled by several fact ors, but, first of all, th e age of the star is important: 
young fast-rotating stars are more act ive (e.g., iTelleschi et ctZ.I . l2005h . The topic of stellar coronal 
loops deserves a review by itself fe.g.. iRosner et all Il985l ) and here only a few relevant issues are 
discussed. A comp l ete an d more recent review of stellar coronae, with an extensive part regarding 
loops, is by iGudell (|200J). In the framework of the solar-stellar connection, it is very important 
the comparison of what we know about the spatially resolved but single solar corona and what 
about the unresolved but numerous stellar coronae, which offer a variety of different environments. 
The lack of spatial resolution inhibits to obtain direct information about the size and appearance 
of the loops, and the general aspect of the corona. We therefore have to rely on indirect evidence. 
One possible approach to get information is to benefit from transient X-ray events, such as flares, 
which prov ide est i mations of t he loop scale length from their dependence on the decay and rise 
timescales (|Realel . l2002al [2003. and Section T4. 1.2 1 for reviews). Detailed hydr odynamic modeling 
can p rovide even more constraints, for instance, on the heat deposition (e.g., iReale et al. . 19881 
20041 ). The study of stellar X-ray flares allo wed, for instance, t o constrain that most stellar flares 

and to infer t he presence both of 
19881 ) and of giant loops 



involve plasma confined in closed structures (|Reale et all 12002 s ). 



loops with size similar to those ob s erved on the Sun (e.g., IReale et 



( Favata et al . 20051 : Getman et al , 2008 ). with length exceeding the stellar radius 

Another approach is to use the entire solar X-ray corona as a template and "Rosetta stone" 
to interpret stellar coronae. A detailed implementation of this approach was devis ed and applied 



extensively using Yohko h data over its entire life, which covers a whole solar cycle (jOrlando et al 



2000; iPeres et a n. I2000I 1. It was shown that the solar corona indeed provides a pattern of com- 
ponents, i.e., quiet structures, active regions, active region cores, flares, which can be identified in 
stellar coronal data and wh i ch can explain ste llar activity giving different weights to the compo- 



nents (jOrlando et al\ . l200ll ; IPeres et al\ . 120041 ) . The method was also applied to describe stellar 



coronae in t erms of loop popul ations and to extract gener al information and constraints on coro- 



nal heating (jPeres et all 120041'). It was applied to flares (jReale et all l200ll ) and to describe the 



evolution of active regions (jOrlando et «Ll , l2004r ). More recently it was shown that a continuous 



unresolved flaring activity may explain the most active coronae, but also that the coronal heating 
appea rs to follow different scaling for quiet regions and for active and flaring regions across the 



cvcle (lArgiroffi et all 120081) . 



Cargill and Klimchuk (|2006l ) realized that the strong hot peaks in the emissio n mea sure- 



temperature distributions in the coronae of some binary stars ( Sanz-Forcada et al. , 2003) 



are 

similar to those expected for an impulsively-heated solar corona. A coronal model comprised of 
many impulsively heated strands shows that the evidence may be compatible with coronae made 
of many very small loops (length under 10 3 km) heated by microflares. 

The recent deeper investigation of solar coronal heating mechanisms through the evidence of 
hot plasma and of variability makes even more important future tighter links t o the study of stellar 
coronae, wh i ch sh ow very strong evidence of very hot steady components (e.g. JSchmitt et al. ■ Il990l : 



Scelsi et aZ.I . [2005V 
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6 Conclusions and Perspectives 



Coronal loops have been the subject of in-depth studies for about 50 years. Since they owe their 
identity to the brightness of the confined plasma, most of the studies have addressed the physics of 
the confined plasma, i.e., its structure, dynamics and evolution. Most of the basic laws that rule 
the confined plasma, such as scaling laws, were developed early after the discovery of loops and 
are now well-established. Although the observational scenario is ever-enriching with the progress 
of the solar coronal missions, a variety of questions remain still open, at all levels, starting from 
loop identification itself. The lack of operative definitions and automatic identification tools, and 
the difficulty to isolate loops from other surrounding and intersecting structures have prevented 
systematic studies on large and unbiased loop populations. Detailed morphological and theoretical 
analyses converge to a fine loop substructuring, which is critical to understand the basic mechanisms 
but is below the resolution limit of present-day instruments. Also the large scale structure of loops 
leaves room for further developments, and in particular the link with the confining magnetic field, 
which is difficult to measure in the corona. One specific question to be addressed is the role of the 
loop tapering in the transition region. 

Looking inside the loops, open questions involve the detailed thermal structure of the con- 
fined plasma. This issue is important to assess the basic loop heating mechanism: a broad 
multi-component thermal distribution would indicate a structured heating, a simpler distribu- 
tion a monolithic mechanism. This question is still open for many reasons; spectroscopic methods 
provide very detailed thermal information, but mostly concentrated in the regime of warm loops. 
Moreover, the methods of analysis do not appear to provide unique answers yet. Also filter ratio 
diagnostics from broad-band X-ray and UV imaging telescopes have not been able to provide con- 
clusive results so far. As a consequence, we are unable, right now, to assess the problem of the 
apparently different nature of warm (TRACE, SoHO/EIT) and hot (Yohkoh/SXT, Hinodc/XRT) 
loops, and whether it is simply a matter of different heating rate or the heating mechanisms are 
radically different. 

The role of the dynamics of plasma confined inside loops is also still under investigation. The 
measurement of plasma motions is made difficult by the possible ambiguity with the apparent 
motion of thermal fronts. It will be important in the future to evaluate the relative weight of the 
different flows, i.e., downflows, evaporation and draining, siphon flows, and their influence on the 
overall loop budget. 

Also the investigation of temporal variations deserves attention. In particular, in narrow band 
instruments thermal variations might be confused with intensity variations, and make the interpre- 
tation difficult. On the other hand, the analysis of emission variations is very important, because 
it can potentially shed light on heating mechanisms based on short impulsive events (nanoflares) 
or on wave- like phenomena (Alfven waves) . 

On the theoretical side, 1-D loop models are well-established and have provided a wealth of 
sound and important results. Today, their evolution consists essentially in their transformation 
into "strand" models, and so a collection of them describes a proper loop. Moreover, loops are 
now seen much more dynamic than they were in the past, either as the site of flows, or of wave 
perturbations, or of heat pulses. So they need more and more time-dependent modeling to ad- 
dress, for instance, the importance of flows, and the relative weight between evaporation-draining 
flows and siphon flows. The modeling is therefore becoming more and more demanding, although 
computing resources are now much more powerful than in the past and although some approaches 
are able to squeeze the spatial dimension, yet maintaining the temporal description. At the same 
time, new data seem to require other model refinements, such as the description of loop expansion 
in the transition region. The improvement of numerical and computing resources is also allow- 
ing more complex and complete modeling of whole loop regions including the 3-D magnetic field 
structure "ab initio" . This approach is very promising and will surely provide important results in 
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the next future, to complement those obtained with basic single loop models. We all look forward 
self-consistent descriptions including the generation of heat from magnetic field rearrangements 
and perturbations. 

Special attention still deserves the problem of what heats coronal loops, which means basically 
what heats the whole corona. This problem has revealed to be particularly difficult, essentially 
because of intrinsic physical reasons, and, in particular, i) the highly effective thermal conduction, 
which inhibits the identification of the heating site, and ii) the expected small scales of the heating 
processes, which require so far prohibitive spatial/temporal resolution. Nevertheless, coronal loops 
remain an obvious excellent laboratory to investigate coronal heating mechanisms, because the 
dense plasma confined therein make them bright and easy to observe. The most recent challenge 
offered by coronal loops is probably the increasing evidence that the thin strands they are made of 
are ignited by small scale, rapid, but intense pulses. Most current efforts are devoted to study this 
aspect both from the theoretical/modeling and from the observational point of view. However, 
alternative explanations are actively explored and strongly proposed lately. 

The study of coronal loops is very alive and is the subject of Coronal Loop Workshops, taking 
place every two years, which are site of debate, inspiration of new investigations, and school for 
young investigators. We look forward further great improvements in our knowledge of coronal 
loops from the new mission Solar Dynamic Observatory. 
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